HD-THEP-96-37 

Equation of state and 
coarse grained free energy 
for matrix models 

J. Berges* 

and 

C. Wetterich^ 



Institut fiir Theoretische Physik 
Universitdt Heidelberg 
Philosophenweg 16 
69120 Heidelberg, Germany 



Abstract 

We investigate phase transitions in three dimensional scalar matrix models, with 
special emphasis on complex 2x2 matrices. The universal equation of state for weak 
first order phase transitions is computed. We also study the coarse grained free 
energy. Its dependence on the coarse graining scale gives a quantitative criterion 
for the validity of the standard treatment of bubble nucleation. 
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1 Introduction 



Matrix models are extensively discussed in statistical physics. Beyond the 0{N) sym- 
metric Heisenberg models (Vector models') they correspond to the simplest scalar field 
theories. There is a wide set of difi^erent applications as the metal insulator transition 
[1] or liquid crystals [2] or strings and random surfaces [3] . . . The universal behavior 
of these models in the vicinity of a second order or weak first order phase transition is 
determined by the symmetries and the field content of the corresponding field theories. 
We will consider here models with U{N) x U{N) symmetry with a scalar field in the 
{N, N) representation, described by an arbitrary complex N x N matrix (p. We do not 
impose nonlinear constraints for ip a priori but rather use a 'classical' potential. This en- 
forces nonlinear constraints in certain limiting cases. Among those, our model describes a 
nonlinear matrix model for unitary matrices or one for singular 2x2 matrices. The uni- 
versal critical behavior does not depend on the details of the classical potential and there 
is no difference between the linear and nonlinear models in the vicinity of the limiting 
cases. We concentrate in this paper on three dimensions, relevant for statistical physics 
and critical phenomena in high temperature field theory. 

The cases N — 2, 3 have a relation to high temperature strong interaction physics. At 
vanishing temperature the four dimensional models can be used for a description of the 
pseudoscalar and scalar mesons for N quark flavors. For = 3 the antihermitean part 
of (fi describes here the pions, kaons, rj and rj' whereas the hermitean part accounts for 
the nonet of scalar 0++ mesons.^ For nonzero temperature T the effects of fluctuations 
with momenta ^ (27rT)^ are described by the corresponding three dimensional models. 
These models account for the long distance physics and obtain after integrating out the 
short distance fluctuations by virtue of dimensional reduction [5, 6, 7]. In particular, the 
three dimensional models embody the essential dynamics in the immediate vicinity of a 
second order or weak first order chiral phase transition [8] . The four dimensional models 
at nonvanishing temperature have also been used for investigations of the temperature 
dependence of meson masses [9, 10]. The simple model investigated in this paper is not 
yet realistic - it neglects the effect of the axial anomaly which reduces the chiral fiavor 
symmetry to SU{N) x SU{N). It also accounts only for meson fluctuations and ignores, 
for example, the temperature dependence of the binding mechanism which gives rise to 
mesons as quark- ant iquark bound states. Nevertheless, it should serve as an interesting 
starting point for a later analysis of more realistic effective three dimensional meson models 
arising in high temperature QCD. For simplicity we will concentrate here on A^ = 2, but 
our methods can be generalized to A^ = 3 and the inclusion of the axial anomaly. 

The case N — 2 also has a relation to the electroweak phase transition in models 
with two Higgs doublets. Our model corresponds here to the critical behavior in a special 
class of left-right symmetric theories in the limit where the gauge couplings are neglected. 
Even though vanishing gauge couplings are not a good approximation for typical realistic 
models one would like to understand this limiting case reliably. 

We are mainly interested in the character and the detailed physics of the phase transi- 

^Scc [4] for a recent phenomenological analysis. 
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tion. Standard (linear) high temperature perturbation theory can give a rehable descrip- 
tion only if the relevant dimensionless couplings remain small near the phase transition. 
This requires two conditions: First, the effective three dimensional couplings evaluated 
for typical momenta ~ (27rT)^ must be small in units of temperature. For the example 
of a quartic coupling this means \{2'kT)/T ^ 1. In high temperature field theory this 
ratio corresponds to the four dimensional coupling, e.g. A(27rT) = A4T. The above con- 
dition is fulfilled whenever the four dimensional (or zero temperature) couplings remain 
perturbative. On the other hand, it never holds for nonlinear matrix models where the 
corresponding ratios diverge. The second condition requires the phase transition to be 
sufficiently strong first order. If we denote by uif. the relevant mass at the critical temper- 
ature, standard perturbation theory turns out to be actually an expansion in X/rric rather 
than \/T. For a given value of A a perturbative expansion can therefore only converge if 
rUc is sufficiently large, rric ^ A4T. In consequence, perturbation theory is not applicable 
for the interesting cases of second order {rric — 0) or weak first order {rric <^ phase 
transitions. 

This property can also be understood in terms of running couplings. One may consider 
the model in presence of an infrared cutoff with typical scale k. Here k may be given by 
external conditions (as a finite volume) or induced by the problem studied (e.g. typical 
external momenta in vertices or a characteristic scale of critical bubbles in problems of 
bubble formation). One may also associate k with a relevant mass away from the critical 
temperature or for nonvanishing sources or external fields. Or else the infrared cutoff can 
be introduced as a purely technical tool. One can define running renormalized couplings 
Ar(A;) by appropriate renormalized n-point functions in presence of the infrared cutoff 
and study the ffow of XR{k) as k is lowered. Often an originally small dimensionless ratio 
A(/c) = \fi{k)/k increases in the course of the evolution. Perturbation theory breaks down 
when \{k) reaches values of order one. One needs nonperturbative methods to follow the 
flow for larger values of A(A;) if k becomes small. This is typically what happens for a 
second order phase transition where one needs at the critical temperature the behavior 
for A; — > 0. In this case the dimensionless coupling \{k) reaches a flxed point A^ which 
is for most simple models substantially larger than one. For the present matrix models 
one wants to know if the phase transition becomes second order in certain regions of 
parameter space. This is equivalent to the question if the system of running couplings 
admits a flxed point which is infrared stable (except one relevant direction corresponding 
to T—Tc). Clearly, this is a nonperturbative question. It can be addressed by an expansion 
in e = 4 — d [11, 12]. However, the reliability of this expansion is not guaranteed a priori 
since for three dimensions e = 1 is not a small parameter. We will employ here a method 
which is based on the general concepts of the Wilsonian renormahzation group [13, 14]. 
More precisely, we study an approximate solution to an exact flow equation for a coarse 
grained free energy. We find that the phase transition for the investigated matrix models 
with N = 2 and symmetry breaking pattern U{2) x f/(2) ^ U{2) is always (fluctuation 
induced) first order, except for a boundary case with enhanced 0(8) symmetry. For a 
large part of parameter space the transition is weak and one finds large renormalized 
dimensionless couplings near the critical temperature. If the running of the couplings 
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towards approximate fixed points (there are no exact fixed points) is sufficiently fast the 
large distance physics looses memory of the details of the short distance or classical action. 
In this case the physics near the phase transition is described by an universal equation of 
state. 

Besides the possible practical applications it is a theoretical challenge to find the 
universal equation of state for weak first order phase transitions. For the second order 
phase transition in 0{N) symmetric vector models (Heisenberg models) this has only 
recently been achieved [15] by the method employed here and confirmed by numerical 
lattice simulations [16]. Whereas for second order transitions the universal equation of 
state can be expressed as a function of only one scaling variable, the additional scale rric 
in a first order transition induces a second dimensionless ratio. The equation of state 
or, equivalently, the effective potential or free energy therefore depends on two different 
scahng variables in the universal region. We have succeeded to solve this problem for the 
present matrix models and present the universal equation of state in sect. 8. 

Another old and challenging problem for first order phase transitions concerns the 
question of the validity of hanger's theory of bubble formation [17]. This requires first 
a meaningful definition of a coarse grained free energy with a coarse graining scale k 
and second the validity of a saddle point approximation for the treatment of fiuctuations 
around the critical bubble. Here only fluctuations with momenta smaller than k must be 
included. We will see that the two issues are closely related. The validity of the saddle 
point approximation typically requires small dimensionless couplings X{k). On the other 
hand we observe for large \{k) that the form of the coarse grained effective potential Uk 
depends strongly on k even for scales k where the location of the minima of Uk is essentially 
independent of k. This means that the lowest order in the saddle point approximation 
(classical contribution) depends strongly on the details of the coarse graining procedure. 
Since the final results as tunneling rates etc. must be independent of the coarse graining 
prescription this is only compatible with a large contribution from the higher orders of 
the saddle point expansion. Section 9 deals with this issue in a quantitative way. 

Our paper is organized as follows. The nonperturbative method we employ relics on 
an exact flow equation for the effective average action which is introduced in section 
2. In section 3 we define the U{2) x U{2) symmetric matrix model and we establish 
the connection to a matrix model for unitary matrices and to one for singular complex 
2x2 matrices. There we also give an interpretation of the model as the coupled system 
of two 5'C/(2)-doublets for the weak interaction Higgs sector. Section 4 is devoted to 
an overview over the phase structure and the coarse grained effective potential Uk for 
the three dimensional theory. The evolution equation for Uk and its scaling form is 
computed in section 5. A method for its numerical solution is discussed in section 6 
which also contains the flow equation for the wave function renormalization constant 
or the anomalous dimension rj. A detailed account on the renormalization group flow 
is presented in section 7. We compute the universal form of the equation of state for 
weak flrst order phase transitions in section 8 and we extract critical exponents and the 
corresponding index relations. The dependence of the coarse grained effective potential 
on the coarse graining scale is studied in detail in section 9. The quantitative analysis is 
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applied to Langer's theory of bubble formation. Section 10 contains the conclusions and 
an outlook. 

2 Effective average action 

The effective average action [18, 19] is a coarse grained free energy with an infrared 
cutoff or a coarse graining scale ~ k. More precisely, for a theory described by a classical 
action S, the effective average action results from the integration of degrees of freedom 
with characteristic momenta larger than k. To be explicit we consider a /c-dependent 
generating functional with N'^ complex scalar fields Xab {a,b— 1 . . . N) 

W,\j] ^IuJdx exp (^-S[x] - A,S[x] + \ j d'x (j:,(x)x"''(x) + xl,{x)3''\x))^ (2.1) 

with S the classical action and jab arbitrary sources. The scale dependence arises from 
the introduction of an additional infrared cutoff term 

^kS[x] = I ^M^)x:b{Q)x''\q). (2.2) 

Without this term Wk becomes the usual generating functional for the connected Green 
functions. Here the infrared cutoff function Rk is required to vanish for /c — > and to 
diverge for A; — > oo and fixed q^. This can be achieved, for example, by the choice 

Mq) = TT-^ (2-3) 

where Z^. denotes an appropriate wave function renormalization constant which will be 
defined below. For fiuctuations with small momenta ^ the cutoff Rj. — Z^k^ 
acts like an additional mass term and prevents their propagation. For ^ k"^ the 
infrared cutoff vanishes such that the functional integration of the high momentum modes 
is not disturbed. The expectation value of x the presence of Afe-S'fx] and j reads 
^ab ^ab y_ 2 SWk[j]/ 5j*i,. We define the effective average action via a Legendre 
transform 

Tfcb] = -Wk[j] + lJd'x {j:,{x)x'\x) + xMf'ix)) - AkSM. (2.4) 

As a consequence, as the scale k is lowered F^ interpolates from the classical action S = 
liuik^oo^k to the standard effective action F = lim^^orfc) i-e. the generating functional 
of IPI Green functions [19]. Lowering k results in a successive inclusion of fluctuations 
with momenta <i and therefore permits to explore the theory on larger and larger 
length scales. One can view F^ as the effective action for averages of fields over a volume 
of size ~ A;""^ and the approach is similar in spirit to the block spin action [13, 14] on the 
lattice. The interpolation property of F^ can be used to 'start' at some high momentum 
scale A where Fa can be taken as the classical or short distance action and to solve the 
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theory by following to /c — > 0. The scale dependence of can be described by an 
exact nonperturbative evolution equation [19, 20, 21, 22] 

|r.M = iTv{(rfM + ii.)-^} (2.6) 

where t = \n(k/A). The evolution is described in terms of the exact inverse propagator 
as given by the second functional derivative of with respect to the fields. The trace 
involves a momentum integration as well as a summation over the internal indices which 
count the 2N'^ real scalar fields contained in cp. The additional cutoff function with a 
form like the one given above renders the momentum integration both infrared (IR) and 
ultraviolet (UV) finite. In particular, the direct implementation of the additional mass 
term Rk — Z^k'^ for <^ k"^ into the inverse average propagator makes the formulation 
suitable for dealing with theories which are plagued by infrared problems in perturbation 
theory. The flow equation (2.5) is compatible with all symmetries of the model and can 
be generalized to include possible local gauge symmetries [23]. The exact renormalization 
group equation can be formulated in many different but formally equivalent ways [14, 24] 
and it may be interpreted as a differential form of the Schwinger-Dyson equations [25] . 

Prom its construction, i.e. the inclusion of fluctuations with characteristic momenta 
larger than a given infrared cutoff ~ A;, the effective average action is the appropriate 
quantity for the study of physics at a scale k. It therefore realizes the concept of a 
coarse grained free energy in the sense of ref. [17]. Such a quantity becomes especially 
desirable for a description of first order phase transitions. In contrast to a second order 
phase transition there is an inherent length scale Iq due to a finite correlation length 
at a first order phase transition. This length scale acts as a physical infrared cutoff. 
The coarse grained effective action with k ~ /q ^ accounts for all fluctuations with 
momenta larger than this physical infrared cutoff and it is the appropriate quantity for 
the study of the physics at the scale /q. To be explicit, one may consider the coarse 
grained effective potential Uk- It is obtained from the coarse grained effective action F^ 
for a constant fleld ip. At a first order phase transition there is a nonzero difference between 
the field expectation value (order parameter) in the symmetric (disordered) and in the 
spontaneously broken (ordered) phase. The two phases correspond to two different minima 
of the coarse grained effective potential Uk and both minima are separated by a potential 
barrier. The coarse grained potential is therefore a nonconvex function whereas the 
standard effective potential U — lim^^o Uk has to be convex by its definition as a Legendre 
transform. The convexity is due to the effect of fluctuations with characteristic length 
scales larger than /q [26]. The study of physical processes such as tunneling or inflation 
usually relies on the nonconvex part of the potential which is discussed in section 9. 

Though the evolution equation (2.5) for the effective average action is exact, it remains 
a complicated functional differential equation. In practice one has to flnd a truncation 
for Tk in order to obtain approximate solutions. An important feature of the exact flow 
equation is therefore its simple and intuitive form which helps to find a nonperturbative 
approximation scheme. The r.h.s. of cq. (2.5) expresses the scale dependence of F^ in 
terms of the exact propagator. Known properties of the propagator can be used as a 
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guide to find an appropriate truncation for tlie effective average action. For a scaiar 
tlieory tfie propagator is a matrix cfiaracterized by mass terms and kinetic terms ~ . 
Tlie mass matrix is given by tfie second derivative of tlie potential C/fc witfi respect to tfie 
fields. In general Z can be a complicated function of the fields and momenta. Wc may 
exploit the fact that the function Z plays the role of a field and momentum dependent 
wave function renormalization. For second order phase transitions and approximately 
for weak first order phase transitions the behavior of Z is governed by the anomalous 
dimension r\. Typically for three and four dimensional scalar theories is small. (In our 
case the relevant value is ^ 0.022 as given by the corresponding index in the 0(8) 
symmetric Hciscnbcrg model). One therefore expects a weak dependence of Z on the 
fields and momenta. We will exploit this in the following. 



3 [/(2) X [/(2) symmetric scalar matrix model 

We consider a C/(2) x [/ (2) symmetric effective action for a scalar field </? which transforms 
in the (2,2) representation with respect to the subgroup SU{2) x SU{2). Here </? is 
represented by a complex 2x2 matrix and the transformations are 

(3.1) 

where U and V are unitary 2x2 matrices corresponding to the two distinct U (2) factors. 

We classify the invariants for the construction of the effective average action by the 
number of derivatives. The lowest order in a systematic derivative expansion [27, 22] of 
Ffc is given by 

Ffe = y d'^x {Uk{^, (^t) + Zkd^^l.d'^^'^'} (a, 6=1,2). (3.2) 

The term with no derivatives defines the scalar potential Uk which is an arbitrary function 
of traces of powers of (p^(p. The most general U{2) x U{2) symmetric scalar potential can 
be expressed as a function of only two independent invariants, 

p = tr((^V) 

1 ^ 2 (3 3) 

T = 2tr((^V--p) =2tr((^V)'-p2. 

Here we have used for later convenience the traceless matrix (p'^cp — ^p to construct the sec- 
ond invariant. Higher invariants, tr (^<^V ~ \p) for n > 2, can be expressed as functions 
of p and T [28]. 

For the derivative part we consider a standard kinetic term with a scale dependent 
wave function renormalization constant Z^. The first correction to the kinetic term would 
include field dependent wave function renormalizations Zk^p.r) plus functions not spec- 
ified in eq. (3.2) which account for a different index structure of invariants with two 
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derivatives. These wave function renormalizations may be defined at zero momentum. 
The next level involves invariants with four derivatives and so on. We define at the 
minimum ''o of Uk and at vanishing momenta q^, 

Zk = Zk{p = Po, T = To; = 0). (3.4) 

The factor appearing in the definition of the infrared cutoff in eq. (2.3) is identified 
with (3.4). The /c-dependence of this function is given by the anomalous dimension 

rj{k) = -^InZfe. (3.5) 

If the ansatz (3.2) is inserted into the flow equation for the effective average action 
(2.5) one obtains flow equations for the effective average potential Uk{p,T) and for the 
wave function renormalization constant Z^. (or equivalently the anomalous dimension 77). 
This is done in sections 5 and 6. These flow equations have to be integrated starting from 
some short distance scale A and one has to specify U^^ and as initial conditions. At 
the scale A, where Fa can be taken as the classical or short distance action, no integration 
of fluctuations has been performed. The short distance potential is taken to be a quartic 
potential which is parametrized by two quartic couplings Aia, A2A and a mass term. We 
start in the spontaneously broken regime where the minimum of the potential occurs at 
a nonvanishing field value and there is a negative mass term at the origin of the potential 

{A > 0), ^ ^ 

^a(p, t) = -filp + 2^iaP^ + 4^2A'r (3.6) 

and Za = 1. The potential is bounded from below provided Aia > and A2A > — 2Aia. 
For A2A > one observes the potential minimum for the configuration ^ah = V^ab 
corresponding to the spontaneous symmetry breaking down to the diagonal U{2) sub- 
group of C/(2) X C/(2). For negative A2A the potential is minimized by the configuration 
ifah = <^SaiSab which corrcspouds to the symmetry breaking pattern U{2) x U(2) — > 
U{1) X f/(l) X f/(l). In the special case A2A = the theory exhibits an enhanced 0(8) 
symmetry. This constitutes the boundary between two phases with different symmetry 
breaking patterns. 

The limits of infinite couplings correspond to nonlinear constraints in the matrix 
model. For Aia 00 with fixed ratio p\/\i\ one finds the constraint ti^ip'^ip) = 2p\/\i\. 
By a convenient choice of Z\ (rescaling of (/?) this can be brought to the form tr{ip''Lp) = 2. 
On the other hand, the limit A2A — > +00 enforces the constraint '^"^Lp = ^trlip'^ip). Com- 
bining the limits Aia 00, A2A 00 the constraint reads = 1 and we deal with 
a matrix model for unitary matrices. (These considerations generalize to arbitrary N.) 
Another interesting limit obtains for Aia = ~|^2A + ^a, > if A2A —00. In this 
case the nonlinear constraint reads (tr (/?''' (/9)^ = tr((/9''"(/9)^ which implies for N — 2 that 
det (f = 0. This is a matrix model for singular complex 2x2 matrices. 

One can also interpret our model as the coupled system of two SU (2)-doublets for the 
weak interaction Higgs sector. This is simply done by decomposing the matrix (pab into 
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two two-component complex fundamental representations of one of the SU (2) subgroups, 
^ah ~^ ^ih) ^2b- The present model corresponds to a particular left-right symmetric model 
with interactions specified by 

We observe that for a typical weak interaction symmetry breaking pattern the expectation 
values of (fi and ip2 should be aligned in the same direction or one of them should vanish. 
In the present model this corresponds to the choice A2A < 0. The phase structure of 
a related model without the term ~ ((/7|<y92)(v'2'f'i) has been investigated previously [29] 
and shows second or first order transitions^. Combining these results with the outcome 
of this work leads already to a detailed qualitative overview over the phase pattern in a 
more general setting with three independent couplings for the quartic invariants {(fi^pi + 
(P2^2)'^, i^i'^i — ¥'29^2)^ and (<^l972)(952<^i)- We also note that the special case A2A — 2Aia 
corresponds to two Heisenberg models interacting only by a term sensitive to the alignment 
between (/?! and Lp2, i.e. a quartic interaction of the form ((/9|(/9i)^+((/92V'2)^+2(v'iV'2)(v'2V^i)- 
The model is now completely specified and it remains to extract the fiow equations for 
Uk and Zk- Before this is carried out in sections 5 and 6 we present an overview over the 
phase structure in the next section. These results are obtained from a numerical solution 
of the evolution equations. 



4 Phase structure 

In this section we consider the U{2) x U (2) symmetric model in three space dimensions. 
The aim is to give an overview of our results concerning the phase structure and the 
effective average potential. We concentrate here on the spontaneous symmetry breaking 
with a residual U{2) symmetry group. This symmetry breaking can be observed for a 
configuration which is proportional to the identity and with (3.3) one finds r = 0. In this 
case we shall use an expansion of Uk{p,T) around r = keeping only the linear term in 
r. This amounts to assuming 

^(p,r = 0) = for n>2. (4.1) 

We will motivate this truncation in section 6 where we present a more detailed analysis. 
We make no expansion of [/fe(p, r) in terms of p. This allows the description of a first 
order phase transition where a second local minimum of Uk{p) = Uk{p,T — 0) appears. 
The p-dependencc also gives information about the equation of state of the system. 

For the considered symmetry breaking pattern the short distance potential Ua given 
in eq. (3.6) is parametrized by positive quartic couplings, 

AiA,A2A>0 (4.2) 

^First order phase transitions and coarse graining have also been discussed in a multi-scalar model 
with Z2 symmetry [30]. 
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and the location of its minimum is given by 

PoA = /^a/Aia. (4.3) 

To study the phase structure of the model we integrate the flow equation for the effective 
average potential C4 (cf. sect. 5, 6) for a variety of initial conditions poA, \a and A2A. In 
particular, for general Aia, A2A > we are able to find a critical value poA = Poc for which 
the system exhibits a first order phase transition. In this case the evolution of Uk leads 
at some scale /c2 < A to the appearance of a second local minimum at the origin of the 
effective average potential and both minima become degenerate in the limit k ^ 0. If 
Po{k) > denotes the /c-dependent outer minimum of the potential (C/^(po) = 0, where 
the prime on Uk denotes the derivative with respect to p at fixed k) at a first order phase 
transition one has 

lim(t/fc(0) - t/fe(po)) = 0. (4.4) 

k—^0 

A measure of the distance from the phase transition is the difference Ska = (poA — Poc)/A. 
If fi\ and therefore poA is interpreted as a function of temperature, the deviation 6k\ is 
proportional to the deviation from the critical temperature Tc, i.e. Sk\ = A{T){Tc — T) 
with A{%) > 0. 

We consider in the following the effective average potential Uk for a nonzero scale 
k. This allows to observe the nonconvex part of the potential (cf. sect. 9). As an ex- 
ample we show in fig. 1 the effective average potential Uk=kf ^oi Aia = Aia/A = 0.1 and 
A2A = A2A/A = 2 as a function of the renormalizcd field ip^ = (pk/2)^/^ with pji = Zk=kfP- 
The scale kf is some characteristic scale below which the location of the minimum po{k) 
becomes essentially independent of k. Its precise definition is given below. We have nor- 
malized Ukf and (fR to powers of the renormalized minimum <por(A;/) = {poR{kf)/2Y^^ 
with poR{kf) = Zk^.po{kf). The potential is shown for various values of deviations from the 
critical temperature or Sk^. For the given examples 6k,a = —0.03, —0.015 the minimum at 
the origin becomes the absolute minimum and the system is in the symmetric (disordered) 
phase. Here (poR denotes the minimum in the metastable ordered phase. In contrast, for 
Ska — 0.04, 0.1 the absolute minimum is located at <Pr/<Por = 1 which characterizes the 
spontaneously broken phase. For large enough 5ka the local minimum at the origin van- 
ishes. For Ska = the two distinct minima are degenerate in height^. As a consequence 
the order parameter makes a discontinuous jump at the phase transition which character- 
izes the transition to be first order. It is instructive to consider some characteristic values 
of the effective average potential. In fig. 2 we consider for Aia = 0.1, A2A = 2 the value 
of the renormalized minimum poR{kf) and the radial mass term as a function of —Ska or 
temperature. In the spontaneously broken phase the renormalized radial mass squared is 
given by (cf. section 5) 

ml{kf) = 2Z^^'poU'^^.ipo), (4.5) 

^ We note that the critical temperature is determined by condition (4.4) in the hmit A; — > 0. Nevertheless 
for the employed nonvanishing scale k = kf the minima of Uk become almost degenerate at the critical 
temperature. 
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Figure 1: The effective average potential Uk=kf as a function of the renormalized field ipR. The potential 
is shown for various values of Ska ^ T,- — T. The parameters for the short distance potential Ua are (1) 
5ka = -0.03, (2) 5ka = -0.015, (3) 6ka = 0, (4) 6ka = 0.04, (5) 5ka = 0.1 and Aia = 0.1, A2A = 2. 

in the symmetric phase the renormahzed mass term reads 

mlR{kf)^Z^lU'^^{Q). (4.6) 

At the critical temperature (^ka = 0) one observes the discontinuity ApoR = PoRikf) 
and the jump in the mass term Auir = mR{kf) — moR^kf) — rrf^ — rnQj^. (Here the 
index 'c' denotes Sk,\ = 0). The ratio Apofi/A is a rough measure for the 'strength' of 
the first order transition. For Apo/j/A ^ 1 the phase transition is weak in the sense that 
typical masses are small compared to A. In consequence, the long-wavelength fluctuations 
play a dominant role and the system exhibits universal behavior, i.e. it becomes largely 
independent of the details at the short distance scale A~^. We will discuss the universal 
behavior in more detail below. 

In order to characterize the strength of the phase transition for arbitrary positive 
values of Aia and A2A we consider lines of constant Apo^/A in the Aia, A2A plane. In flg. 
3 this is done for the logarithms of these quantities. For fixed A2A one observes that the 
discontinuity at the phase transition weakens with increased Aia- On the other hand for 
given AiA one finds a larger jump in the order parameter for increased A2a- This is true 
up to a saturation point where Apqr/A becomes independent of A2a- In the plot this can 
be observed from the vertical part of the line of constant ln(Apoij/A). This phenomenon 
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Figure 2: The minimum pofl and the mass term ma in units of the momentum scale A as a function 
of —Ska or temperature (Aia = 0.1, A2A = 2, k = kf). For Ska = one observes the jump in the 
renormalized order parameter Apofl and mass AmR. 

occurs for arbitrary nonvanishing Apon/A in the strong A2A coupling limit and is discussed 
in section 7. 

In the following we give a detailed quantitative description of the first order phase 
transitions and a separation in weak and strong transitions. We consider some charac- 
teristic quantities for the effective average potential in dependence on the short distance 
parameters Aia and A2A for 6ka = 0. We consider the discontinuity in the renormalized 
order parameter ApoR and the inverse correlation lengths (mass terms) and ttIq^ in 
the ordered and the disordered phase respectively. Fig. 4 shows the logarithm of Apoii in 
units of A as a function of the logarithm of the initial coupling A2a- We have connected 
the calculated values obtained for various fixed Aia = 0.1, 2 and Aia = 4 by straight lines. 
The values are listed in table 1. For A2A/A1A ^ 1 the curves show constant positive slope. 
In this range Apoi? follows a power law behavior 

Apoi^ = i?(A2A)^ ^ = 1.93. (4.7) 

The critical exponent 9 is obtained from the slope of the curve in fig. 4 for A2A/A1A <S 1. 
The exponent is universal and, therefore, does not depend on the specific value for Aia. 
On the other hand, the amplitude R grows with decreasing Aia- For vanishing A2A the 
order parameter changes continuously at the transition point and one observes a second 
order phase transition as expected for the 0(8) symmetric vector model. As A2A/A1A 
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Figure 3: Lines of constant jump of the renormalized order parameter Apoi?/A at the phase transition 
in the In(AiA), ln(A2A) plane. The curves correspond to (1) \n{ApoR/A) = —4.0, (2) ln(Apoii/A) = —4.4, 
(3) ln(Apoij/A) = -10.2, (4) ln(Apoij/A) = -14.3. 

becomes larger than one the curves deviate substantially from the linear behavior. The 
deviation depends on the specific choice of the short distance potential. For A2a/ Aia ^ 1 
the curves flatten. In this range Apoi? becomes insensitive to a variation of the quartic 
couphng A2A- 

In addition to the jump in the order parameter we present the mass terms and 
ttIq^ which we normalize to ApoR- In fig- 5 these ratios are plotted versus the logarithm 
of the ratio of the initial quartic couplings A2a/Aia- Again values obtained for fixed 
AiA = 0.1, 2 and Aia = 4 are connected by straight lines. The universal range is set 
by the condition m^/ApoR — const (equivalently for mQj^/Apoji). The universal ratios 
are m^/ApoR = 1-69 and mgjij/ApoK = 1-26 as can be seen from table 1. For the given 
curves universality holds approximately for A2a/Aia 1/2 and becomes 'exact' in the limit 
A2A/A1A 0. In this range we obtain 



The universal features of the system are not restricted to the weak coupling region of 
A2A- This is demonstrated in fig. 5 for values up to A2A — 2. The ratios m^/Apoi? and 
mQ^/ApoR deviate from the universal values as A2A/ Aia is increased. For fixed A2A a larger 
AiA results in a weaker transition concerning Apo^/A. The ratio m^/Apo_R increases with 
AiA for small fixed A2A whereas in the asymptotic region, A2A/A1A ^ 1, one observes from 



m% = S{X2Ay, 



m^^ = ^(A2a)'. 



(4.8) 
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Figure 4: The logarithm of the discontinuity of the rcnormahzed order parameter Apofl/A as a function 
of ln(A2A)- Data points for fixed (1) Aia = 0.1, (2) Aia = 2, (3) Aia = 4 are connected by straight hnes. 

fig. 5 that this tendency is reversed and mfj/Apofl, f^oRl ^Pm start to decrease at about 
AiA ~ 2. 

In summary, the above results show that though the short distance potential ?7a indi- 
cates a second order phase transition, the transition becomes first order once fiuctuations 
are taken into account. This fluctuation induced flrst order phase transition is known in 
four dimensions as the Coleman- Weinberg phenomenon [31] . Previous studies of the three 
dimensional U(2) x U{2) symmetric model using the e-expansion [32] already indicated 
that the phase transition should be a fluctuation induced flrst order transition. The valid- 
ity of the e-cxpansion for weak flrst order transitions is, however, not clear a priori since 
the expansion parameter is not small - there are known cases where it fails to predict 
correctly the order of the transition [33] . The question of the order of the phase transition 
has been addressed also in lattice studies [34] and in high-temperature expansion [35]. All 
studies are consistent with the flrst order nature of the transition and with the absence 
of nonperturbative infrared stable flxed points. It is, however, notoriously difficult to 
distinguish by these methods between weak first order and second order transitions. Our 
method gives here a clear and unambiguous answer and allows a detailed quantitative 
description of the phase transition. The universal form of the equation of state for weak 
first order phase transitions is presented in section 8. 

In the following we specify the scale for which we have given the effective average 
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mm 


A 






A 






0.1 


0.005 


0.386 X 10 


1.69 


1.26 


2 


2 


0.392 X 10 ^ 


1.66 


1.24 


0.1 


0.01 


0.158 X 10 ^ 


1.68 


1.26 


2 


6 


0.230 X 10 ^ 


1.68 


1.25 


0.1 


0.05 


0.461 X 10 * 


1.66 


1.23 


2 


15 


0.505 X 10 ^ 


1.80 


1.35 


0.1 


0.1 


0.249 X 10 


1.58 


1.17 


2 


30 


0.649 X 10 ^ 


1.91 


1.45 


0.1 


0.2 


0.193 X 10 


1.34 


0.992 


2 


70 


0.712 X 10 ^ 


2.01 


1.60 


0.1 


0.5 


0.316 X 10 ^ 


0.772 


0.571 


2 


150 


0.699 X 10 ^ 


2.02 


1.65 


0.1 


1 


0.145 


0.527 


0.395 


2 


350 


0.685 X 10 ^ 


2.03 


1.68 


0.1 


2 


0.341 


0.455 


0.360 


4 


0.2 


0.298 X 10 * 


1.69 


1.26 


0.1 


5 


0.547 


0.450 


0.414 


4 


2 


0.213 X 10 


1.70 


1.27 


0.1 


10 


0.610 


0.462 


0.490 


4 


12 


0.195 X 10 ^ 


1.80 


1.35 


2 


0.02 


0.523 X 10 


1.69 


1.26 


4 


30 


0.302 X 10 


1.89 


1.43 


2 


0.1 


0.121 X 10"^ 


1.69 


1.26 


4 


70 


0.355 X 10-1 


1.96 


1.49 


2 


0.3 


0.101 X IQ-^ 


1.69 


1.25 


4 


150 


0.369 X 10-1 


1.98 


1.55 


2 


1 


0.104 X 10-2 


1.66 


1.25 


4 


350 


0.372 X 10-1 


1.97 


1.57 



Table 1 : The discontinuity in the renormaUzed order parameter ApoR and the critical inverse correlation 
lengths m^j and mQ^j in the ordered and the disordered phase respectively. For small A2A/A1A the ratios 
m^/Apofl and mg^j/Apofl become universal. 



potential Uk- We observe that Uk depends strongly on the infrared cutoff k as long as k 
is larger than the scale k2 where the second minimum of the potential appears. Below 
k2 the two minima start to become almost degenerate for T near and the running of 
Po{k) stops rather soon. The nonvanishing value of k^ induces a physical infrared cutoff 
and represents a characteristic scale for the first order phase transition. We stop the 
integration of the flow equation for the effective average potential at a scale kf < k2 
which is determined in terms of the curvature (mass term) at the top of the potential 
barrier that separates the two local minima of Uk at the origin and at po{k). The top of 
the potential barrier at psik) is determined by 

U'M = (4.9) 

for < Pb(/c) < Po{k) and for the renormalized mass term at psik) one obtains 

ml,R{k) = 2Z^'pBU'apB) < 0. (4.10) 
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We fix our final value for the running by 

kj - \ml,Rikf)\ 



0.01 (4.11) 



For this choice the coarse grained effective potential Ukf essentially includes all fluctua- 
tions with momenta larger than the mass |mB^ij| at the top of the potential barrier. It is a 
nonconvex function which is the appropriate quantity for the study of physical processes 
such as tunneling or inflation. The nonconvex part of Uk is considered in detail in section 
9. There we will also discuss that the appropriate choice for a coarse graining scale k is 
often far from obvious. 



5 Scale dependence of the effective average potential 

In this section we derive the flow equation for the effective average potential Uk- We point 
out a related investigation of the four dimensional SU {N) x SU (N) symmetric linear sigma 
model coupled to fermions that has been studied previously within the framework of the 
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effective average action [28]. Tliere tlie flow equations for a polynomial approximation of 
Uk can be found for arbitrary dimension d. To study the equation of state of the three 
dimensional theory we keep here the most general form for Uk- 

We define Uk by evaluating the average action for a constant field with = QU^ 
where Q denotes the volume. To evaluate the r.h.s. of (2.5) with the ansatz (3.2) we 
expand around a constant background configuration. With the help of the U{2)xU (2) 
transformations the matrix field </? can be turned into a standard diagonal form with 
real nonnegative eigenvalues. Without loss of generality the evolution equation for the 
effective potential can therefore be obtained by calculating the trace in (2.5) for small field 
fluctuations Xab around a constant background configuration which is real and diagonal, 

Vah = VJah , V*a = Va- (5-1) 

We separate the fluctuation fleld into its real and imaginary part, Xab = -^{XRah + 
iXiah) and perform the second functional derivatives of with respect to the eight real 
components. For the constant configuration (5.1) it turns out that F^^^ has a block 
diagonal form because mixed derivatives with respect to real and imaginary parts of 
the field vanish. The remaining submatriccs Sx'r^x'r ^'^'^ ^^^kl ^x'f^X'f can be 

diagonalized in order to find the inverse of F^, + Rk under the trace occuring in eq. (2.5). 
Here the momentum independent part of F^ defines the mass matrix by the second 
functional derivatives oiUk- The eight eigenvalues of the mass matrix are 

{Mtf = C/^ + 2(p±(p^-r)V2)9.C/,, 

{M^f = U'k±2T'/'drUk 

corresponding to second derivatives with respect to xi and 

(Mi)' = (Mf)^ 

(M^)2 = U'k + pU',! + 2pdrUk+ATdrU'k + 4pTd^,Uk 

± {r {Ul! + AdrUk + ipdrU'k + Ard^M^ ^^"^^ 

+ (p2 - r) (C/^ - 2d,Uk - Ard'^UkfY^' 

corresponding to second derivatives with respect to xr- Here the eigenvalues are expressed 
in terms of the invariants p and r using 

<pl-l{p + r'^% <pl = lip-r'/') (5.4) 

and we adopt the convention that a prime on Uk{p, r) denotes the derivative with respect 
to p at fixed r and k and d'!^'Uk = c?"?7fe/(5r)'*. 

The fiow equation for the effective average potential is simply expressed in terms of 



(5.2) 
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the mass eigenvalues 

d If d'^Q 



+ , -.^^2 + , .L+.„^^^2 (5-5) 



{Pk{q) + {M+{p,T)y Pk{q) + {M,-{p,T)y P^iq) + {M} {p, T)y 
111 



Pfc(g) + (M2-(p,r))2 P,(g) + (M4+(p,r))2 P,(g) + {M^{p,T)yj- 

On the right hand side of the evolution equation appears the (massless) inverse average 
propagator 

Pfe(g) = Zkq' + Rk{q) = ^ _^'_^,2/fe2 (5.6) 

which incorporates the infrared cutoff function given by eq. (2.3). The only approxima- 
tion so far is due to the derivative expansion (3.2) of Tk which enters into the flow equation 
(5.5) through the form of P^.. The mass eigenvalues (5.2) and (5.3) appearing in the above 
flow equation are exact since we have kept for the potential the most general form Uk{p, r). 

Spontaneous symmetry breaking and mass spectra 

In the following we consider spontaneous symmetry breaking patterns and the corre- 
sponding mass spectra for a few special cases. For the origin at ipah = all eigenvalues 
equal [/^(O, 0). If the origin is the absolute minimum of the potential we are in the 
symmetric regime where all excitations have mass squared [/^(O, 0). 

Spontaneous symmetry breaking to the diagonal U{2) subgroup of U{2) x U{2) can 
be observed for a field configuration which is proportional to the identity matrix, i.e. 
^ab = fSah- The invariants (3.3) take on values p = 2ip^ and r = 0. The relevant 
information for this symmetry breaking pattern is contained in Uk{p) = Uk{p,T = 0). In 
case of spontaneous symmetry breaking there is a nonvanishing value for the minimum 
Po of the potential. With [/^(po) — one finds the expected four massless Goldstone 
bosons with {M^Y — {M^Y — (-^;r)^ — 0- addition there are three massive scalars in 
the adjoint representation of the unbroken diagonal SU{2) with mass squared {M^Y = 
(M3+)2 = (M4")2 = ApodrUk and one singlet with mass squared (M/)^ = 2pof/fc- The 
situation corresponds to chiral symmetry breaking in two fiavor QCD in absence of quark 
masses and the chiral anomaly. The Goldstone modes are the pseudoscalar pions and the 
T] (or rj') , the scalar triplet has the quantum numbers of qq and the singlet is the so-called 
(7-field. 

Another interesting case is the spontaneous symmetry breaking down to a residual 
U{1) X U{1) X U{1) subgroup of U{2) x U{2) which can be observed for the configuration 
^ab — 'P^ai^ab (p — , T — (fi^ — p^) . Corresponding to the number of broken generators 
one observes the five massless Goldstone bosons (Mj^)^ = (M2")^ = (-^3^)^ = for the 
minimum of the potential at U^. + 2podrUk = 0. In addition there are two scalars with 
mass squared (M^)^ = {M^f = K - '^PodrUk and one with {M^f = U'^ + 2poU'^ + 
6podrUk + 8pldrUi + 8pld^,Uk. 
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We finally point out the special case where the potential is independent of the second 
invariant r. In this case there is an enhanced 0(8) symmetry instead of U{2)xU (2). With 
9"C/fe = and Ul{po) ~ one observes the expected seven massless Goldstone bosons and 
one massive mode with mass squared 2pqU'^. 

Scaling form of the flow equation 

For the 0(8) symmetric model in the hmit A2A = one expects a region of the 
parameter space which is characterized by renormalized masses much smaller than the 
ultraviolet cutoff or inverse microscopic length scale of the theory. In particular, in the 
absence of a mass scale one expects a scaling behavior of the effective average potential 
Uk- The behavior oiU^ at or near a second order phase transition is most conveniently 
studied using the scaling form of the evolution equation. This form is also appropriate for 
an investigation that has to deal with weak first order phase transitions as encountered 
in the present model for A2A > 0. The remaining part of this section is devoted to the 
derivation of the scaling form (5.15) of the flow equation (5.5). 

In the present form of eq. (5.5) the r.h.s. shows an explicit dependence on the scale k 
once the momentum integration is performed. By a proper choice of variables we cast the 
evolution equation into a form where the scale no longer appears explicitly. We introduce a 
dimensionless potential Uk — k~'^Uk and express it in terms of dimensionless renormalized 
fields 

p = Zkk^-^p, 



Zlk^-^^T . 



(5.7) 



The derivatives of arc given by 

9>f ^ {p, f) = Z^-2n-m^(2n+m-l)d-4n-2mQn^M (^^ ^) _ (53) 

(Note that m^™"* denotes m derivatives with respect to p at fixed f and k, while uj^"^^ 
denotes m derivatives with respect to p at fixed r and k). With 

d 

^Mfc(p,f)|p,f = -duk(p,f) + (d-2 + r])pu'k{p,f) + (2d-A + 2r])fdfUk(p,f) 

Q (5-9) 

+k-'-U,{p{p),T{f))^^^^ 

one obtains from (5.5) the evolution equation for the dimensionless potential. Here the 
anomalous dimension 77 arises from the t-derivative acting on and is given by eq. (3.5). 
It is convenient to introduce dimensionless integrals by 



2 J {2ny dt Pk + Zkk^(jj 
where 



I'^ico) - rj^{u) (5.10) 



vt = ( ^ ) . (5.11) 
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The explicit form of /q and Iq reads 



/oV) = -l^'dyy'^+'^[y{l + r{y))+u;] 

^ 1 /'°° d 

^oM = o/ o??/y^r(y)[y(l + r(y))+a;]-' 
2 JO 



(5.12) 



2 Jo 

with the dimensionless infrared cutoff function 



The behavior of these functions is discussed in section 6 where we also consider their 
derivatives with respect to cu. Using the notation 

/oV;^) = ^oH-^/oH (5.14) 

one obtains the flow equation for the dimensionless potential 

d ^ ^ 

w:Uk{p, t) = -dukip, T) + {d-2 + 7])pu',^{p, T) + {2d-A + 2ri)TdfUk{p, r) 
at 

+Avdt {{mtip, f)f; v) + ((mr(p, f))^; 77) + 2vX ((m+(p, f))^; 77) (5-15) 
+2v4 {{m^{p, f)f; ri) + 2vX ( f))^; 77) + 2v4 {{m^ip, f)f; v) 
where the dimensionless mass terms are related to (5.3) according to 

{mfip, t))2 = ^ ^. (5.16) 

Eq. (5.15) is the scaling form of the flow equation we are looking for. For a f-indepcndent 
potential it reduces to the evolution equation for the 0(8) symmetric model [27, 22]. The 
potential Uk at a second order phase transition is given by a /^-independent (scaling) solu- 
tion duk/dt — [27, 22]. For this solution all the /c-dependent functions on the r.h.s. of eq. 
(5.15) become independent of k. For a weak first order phase transition these functions 
will show a weak A';-dcpcndcnce for k larger than the inherent mass scale of the system 
(cf. section 7). There is no particular advantage of the scaling form of the flow equation 
for strong first order phase transitions. 



6 Solving the flow equation 

Eq. (5.15) describes the scale dependence of the effective average potential by a non- 
linear partial differential equation for the three variables t, p and f. A major difficulty for 
its analytical study is that the integral /Q(a;;?7) (cf. eq. (5.12)) can be done analytically 
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only for certain limits of the arguments. The complicated form of the equation therefore 
suggests a numerical solution. We will use a method that relies on a simultaneous ex- 
pansion of the potential around a number of field values Pi, fj, i — 1, . . . ,1, j — 1, . . . ,1' 
for given numbers /, The expansions around different points are matched to obtain 
the general field dependence of the potential. As a consequence we cast the partial dif- 
ferential equation (5.15) into a system of ordinary differential equations. The method is 
developed in [37] for a computation of the critical equation of state for 0{N) symmetric 
models [15]. The generahzation to the present model is described below. This approach 
has some favorable aspects. The main advantage is that it allows the integration of the 
differential equations using a standard Runge-Kutta algorithm without the occurance of 
numerical instabilities. (See e.g. [37] for instability problems with standard discretised 
versions of partial differential equations in the context of flow equations) . The coupled set 
of ordinary differential equations describes the flow of couplings deflned as derivatives of 
the potential at given points, e.g. at the minimum of the potential. These couplings often 
allow direct physical interpretation and some of their properties can be read off from the 
analytic structure of their flow equations. We will exploit this fact to explain the results 
we obtain from the numerical solution in section 7. 

We concentrate in the following on spontaneous symmetry breaking with a residual 
U (2) symmetry group. As we have already pointed out in section 3 this symmetry breaking 
can be observed for a configuration which is proportional to the identity and we have f = 0. 
In this case the eigenvalues (5.2) and (5.3) of the mass matrix with (5.16) are given by 

{mty = (mg )2 = u'k , 

{mty = (mj)2 = < + ApdfUk , (6.1) 
< + '2pK 

and on the r.h.s. of the partial differential equation (5.15) for Uk{p) = Uk{p,f = 0) only 
the functions u'^ip), u'l{p) and dfUk{p) appear. At fixed p — pi the fc-dependence of Uk 
is then determined by the couplings u'i.{pi), u'l{pi) and dfUkipi). Wc determine these 
couplings through the use of flow equations which are obtained by taking the derivative 
in eq. (5.15) with respect to p and f evaluated at p = pj, f = 0. These flow equations for 
u'l^, u'l and dfUk involve also higher derivatives of the potential as m'^', u^^\ dfu'y. and dfu'^. 
The procedure will be to evaluate the flow equations for u'^ and u'l. at different points 
Pi, f = for i = 1, . . . , Z and to estimate the higher p-derivatives appearing on the right 
hand side of the flow equations by imposing matching conditions. The same procedure 
can be applied to dfUj. and in order to obtain an equivalent matching we also consider 
the flow equation for dfu'^,. One could proceed in a similar way for d'^Uk which appears 
on the right hand side of the evolution equation of dfUk and so on. However, since we 
are interested in the p-dependence of the potential at f = we shall use a truncated 
expansion^ in f with 

a>fc(p, f - 0) = for n > 2. (6.2) 

^In principle one could also consider points with f 7^ in the neighborhood of f = and use the 
additional information to estimate the higher f-derivatives as it is done for the higher p-derivatives. 



(mi)2 = 
{mtf = 
{mtf = 
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In three space dimensions the neglected (p-dependent) operators have negative canonical 
mass dimension. We make no expansion in terms of p since the general p-dependence 
allows a description of a first order phase transition where a second local minimum of 
Uk{p) appears. The approximation (6.2) only affects the flow equations for drU^ and dfu'^. 
(cf. eqs. (6.9) and (6.10)). The form of the flow equations for u'f, and u'l, is not affected by 
the truncation (cf. eqs. (6.7) and (6.8)). From m'^ we obtain the effective average potential 
Uk by simple integration. We have tested the sensitivity of our results for u'^, to a change 
in dfUk by neglecting the p-dependence of the f-derivative. We observed no qualitative 
change of the results. We expect that the main truncation error is due to the derivative 
expansion (3.2) for the effective average action. 
To simplify notation we introduce 



e(p) = <(P,7^ = 0), 
Ai(p) = <(p,f = 0), 
A2(p) = 49fitit(p,f = 0). 



(6.3) 



Higher derivatives are denoted by primes on the p-dependent quartic 'couplings', i.e. 
\\ — u'l', A2 = dfu'k etc. It is convenient to introduce functions l'^{uj]r]) that can be 
related to [^{uj^ri) (5.12) by differentiation with respect to the mass argument: 



li{uj-r)) = Ifiuj) - riliiu;) 
1 d 



(6.4) 



for n > 1. 



The explicit form of and /^^ is given in eq. (A.l) in the appendix. For uj > —1 they 
are positive, monotonically decreasing functions of uj. In leading order and vanish 
~ for arguments a; » 1. They introduce a 'threshold' behavior that accounts 

for the decoupling of modes with mass squared larger than the infrared cutoff Z^A;^. For 
vanishing argument they are of order unity. As u —1 the functions exhibit a 

pole for d < 2(n + 1). The pole structure is discussed in the appendix. We also use two- 
parameter functions ^^(a;!, a;2; r/) [28]. For rii = n2 = I their relation to the functions 
l'^{LO]r]) can be expressed as 



1 



UJ2 - UJi 

li{uj-ri) 



77) - l'i{u2] rf) for ui ^ LO2, 



(6.5) 



and 



C+l,n2(^l'^2;?7) 



1 d 



C,n2K>'^2;r?), li^^^Ai^uUJ2]V) = ltnS^2,^l\ri). (6.6) 
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With the help of these functions the scale dependence of e is described by 



de 
di 



(-2 + r))e+(d-2 + ri)p\^ - ^Vd{\i + A2 + PKYHe + 9X2] v) 
-2vd{3X, + 2pX'X{e + 2pAi; r]) - Sv^Xiltie; rj) 



(6.7) 



and for Ai one finds 
dXi 



dt 



= (d - 4 + 27y)Ai + {d-2 + 7])pX[ 



+6va [(Ai + A2 + pX'^f4{e + PA2; v) - (A'l + 2A'2 + pX^e + pX^; r/) 
+2vd [(3Ai + 2pX\)Hi{e + 2pAi; 1]) - (5A; + 2pX';)lf{e + 2pAi; tj) 
+8vd [{X,)Y^{e;rj)-X[li{e;rj) 

Similarly the scale dependence of A2 is given by 

dX2 



(6.8) 



dt 



(d - 4 + 2r])X2 + (d-2 + ri)pX:^ - Avd{X2fl{^{e + pA2, e; r)) 

+2vd [3(A2)' + I2A1A2 + 8pA'2(Ai + A2) + Ap\X'^f] /(^(e + pA2, e + 2pAi; v) 
-UvdKliie + PA2; ry) - 2T;d(5A'2 + 2pX^)li{e + 2pAi; 77) 
+2vd [(A2)%'(e;77)-4A'2Zf(e;?7)' 



(6.9) 



and for A2 it reads 



9A^ 
9i 



{2d-Q + ?>7])X'^ + {d-2 + r])p>:J^ 

+^Vd{X2f [(Ai + A2 + pX'^)lii{6 + pX2, e; 77) + XJ^^i^ + pX2, e; 77) 
-2vd[3A2(4Ai + A2) + 4pA^(2Ai + 2A2 + PA2)] [(^i + A2 + pX'^) 
^2,i(e + pA2, e + 2pXi; r,) + (3Ai + 2pX[)ll^{e + pA2, e + 2pAi; ^) 
+4^;d[A'2(7A2 + lOAi) + 6A2A'i + 4p(A'2'(Ai + A2 + PA2) + A2(2A'2 + A'^)) 
^1,1 (e + pA2, e + 2pAi; v) - 8vdX2X!^li^^{e + pA2, e; 77) 
+2^;rf(3Ai + 2pAi)(5A'2 + 2pX'i)li{e + 2pAi; r)) + 14T;rfA'2(Ai + A2 + PA2) 
li{e + pAs; f]) - 2vd{7X'^ + 2pA'2")^f (e + 2pAi; ^7) - 14^;rfA'2'/f (e + pA2; v) 
-Avd{X2fXili{e; r]) + 4t;,A'2(2Ai + A2)/^(e; r?) - SvdX'^lfie; r]). 

(6.10) 

Wc evaluate the above flow equations at different points pi for i = 1, . . . , / and use a 
set of matching conditions that has been proposed in ref. [37]. The generalization of 
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(6.11) 



these conditions to the present model is obtained by considering fourth order polynomial 
expansions of e(p) and X2{p) around some arbitrary point p,, 

{X2UP) = A2, + A'2,(p-p.) + ^A'2',(p-p,f + 1a^';,(p-pO' 
with ej — e{pi), A2,i = X2{pi) etc. Using similar expressions for 

(Ai).(p) = §z{eUp) , (X'Mp) = Ypi^^Up) (6-12) 

the matching is done by imposing continuity at half distance between neighboring expan- 
sion points, 

W,(^)=(e)„(^i^) , (AO.(^) = (AO„(^ 

Pi + Pi+1 



(A.).(^) = (A,)„(^) , (Aa.(^) = (Aa 



j+1 



2 

(6.13) 

for i = 1 — 1 and 

(Kh (^^) = (A'.).« (^^) . (A3, (^^) = (A3,« (^^ 

(6.14) 

for the initial and end points, J = 1 and j = I — 1. Together these 4(/ — 1) conditions 
(6.13) for all / — 1 > 2 intermediate points and the four conditions (6.14) make up two 
independent algebraic systems of each 21 equations. From the first set of equations one 
obtains a unique solution for A'^ ^ and A^ j. The second set is identical in structure and A2 j, 
X2 .1 can be obtained from the solutions for A'^ ^ and A'/ ^ with the substitutions ej — > A2J, 
Aij — > A2J, X'l j — > X2J, Xi j A2J for j = 1, . . . ,/. With the help of these algebraic 
solutions we eliminate X[{p), A'/(p), A2(p) and X'/^i^p) in the flow equations (6.7) - (6.10) 
for all / points pi.^ Therefore, equations (6.7) - (6.10) are turned into a closed system of 

^The algebraic solutions A'^ ^ and A" ^ (equivalently Aj j and \2\) do incorporate information from the 
whole range of points pj with j = 1, . . . ,1. It is a feature of the matching conditions (6.13), (6.14) that 
the contributions from points Pj^i to \'i j, A" j rapidly decrease with increasing |« — i|. For equal spacings 
between neighboring expansion points contributions from points pj with j > i + 1 (j < i — 1) are typically 
suppressed by a factor ;^10~l'~-'l+^ as compared to the contribution from the nearest neighbor point 
(/5i_i). As a consequence solutions X'^ j, A" ^ for inner points with 1 <^ i <^ I become independent from 
boundary points. We observe approximate translational invariance for inner point solutions, i.e. X'-^ 
and A"^_|_„ are approximately obtained from the solutions A'^ ^ and A" ^ with the substitutions Cj — » ej±„, 
Ai J — > Xi j±n for 1 < J < / if z and i ± n are sufficiently far away from the boundaries. The decoupling 
from distant points and the translational invariance for inner points can be used to obtain approximate 
expressions which become useful if a large number of expansion points is considered. We use the exact 
algebraic solution for I = 10 points. For ^ > 10 we apply the approximate translational invariance to 
generate from A'j^ 5, X'^ q additional solutions A^ 5^211 ''^'1,6+21 for i — 1, ...,(/ — 10)/2 with I even and 
equivalently for A" 5, A" g. With A'^ ;_3, . . . , A'^ ; and A'/;_3, . . . , A"j from the calculation with 10 points 
one obtains the desired generalization. We have used runs with different choices of I in order to check 
the stability of the numerical results. 
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4/ ordinary differential equations for tlie unknowns e(pi), \i{pi), ^2{Pi) and X2{pi)- 

If there is a minimum of the potential at nonvanishing k = powe use expansion points 
that are proportional to the minimum, i.e. pi — with i — 1, . . . ,1 and fixed integer 
n. The condition e(«;) = can be used to obtain the scale dependence of K{k): 



P=K 



= -{d-2 + v)K + 6v, (l + ^^^;^^^^) It {^H^y.ri) (6.15) 

(2k\' (k) \ 
3 + -J^j (2«Ai(«;); r?) + 8t;,/f (0; 7]). 

To make contact with /3-functions for the couplings at the potential minimum k we point 
out the relation 

^^ = ^Ip^. + AU \^)^^- (6.16) 
Similar relations hold for e(pj), Ai(pj) etc., e.g. 

de(pj) Se, i — 1 ,_,dfi; 

dAi(pj) _ dXi , i~l\//~Nd« 
di " 9t '^=^^+ n -^^^^^^dr 

We integrate the 4Z — 1 difTerential equations (6.7) - (6.10) for the couplings e(pi), Ai(pi), 
\2{Pi) and A2(pi) (with d/dt replaced by d/dt according to (6.17)) and the one for k (6.15) 
with a fifth-order Runge-Kutta algorithm using the embedded fourth-order method for 
precision control. The general p-dcpcndcnce is recovered by patching the simultaneous 
expansions around different points at half distance between neighboring expansion points. 

It remains to compute the anomalous dimension r] defined in (3.5) which describes the 

scale dependence of the wave function renormalization Zk- Wc consider a space dependent 
distortion of the constant background field configuration (5.1) of the form 

<Pah{x) = <Pa5ah + [^ife'''^'' + 5ip*e'^''] E,ft. (6.18) 
Insertion of the above configuration into the parametrization (3.2) of yields 

Zk = Zk{p,T,Q'^ = 0) = ^ / lim |jy=o. (6.19) 

To obtain the flow equation of the wave function renormalization one expands the effective 
average action around a configuration of the form (6.18) and evaluates the r.h.s. of eq. 
(2.5). This computation has been done in ref. [28] for a 'Goldstone' configuration with 



-'ab 



8al8h2 - 8a2hl (6.20) 
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and (faSab = ^Sab corresponding to a symmetry breaking pattern with residual U{2) 
symmetry. The result of ref. [28] can be easily generalized to arbitrary fixed field values 
of p and we find 

r]{k) = 4^p[4(Ai)^m^,2(e,e + 2pAi;^) + (A2)'m^,2(e,e + pA2;r?)]. (6.21) 
The explicit form of the 'threshold' function 

777,2^2 ("^l' '^2;^) = m2^2i'^l, (JJ2) - r]'!^^^^^'^^, UJ2) (6.22) 

can be found in refs. [27, 28]. For vanishing arguments the functions 777.22 ^22 
of order unity. They are symmetric with respect to their arguments and in leading order 
77722(0,60') ~ 77722(0,0^) ~ o;"^ for a; 1. According to eq. (3.4) we use p = n to define 
the uniform wave function renormalization 

Zk = Zk{K). (6.23) 

We point out that according to our truncation of the effective average action with eq. (6.21) 
the anomalous dimension 77 is exactly zero at p = 0. This is an artefact of the truncation 
and we expect the symmetric phase to be more affected by truncation errors than the 
spontaneously broken phase. We typically observe small values for 7; (A;) = — d(lnZfc)/dt 
(of the order of a few per cent). The smallness of 77 is crucial for our approximation 
of a uniform wave function renormalization to give quantitatively reliable results for the 
equation of state. For the universal equation of state given in sect. 8 one has 77 = 0.022 
as given by the corresponding index of the 0(8) symmetric 'vector' model. 



7 Renormalization group flow 

To understand the detailed picture of the phase structure presented in section 4 we will 
consider the fiow of some characteristic quantities for the effective average potential as the 
infrared cutoff k is lowered. We will always consider in this section the trajectories for the 
critical 'temperature', i.e. = 0, and we follow the flow for different values of the short 
distance parameters Aia and A2A. The discussion for sufficiently small is analogous. 
In particular, we compare the renormalization group flow of these quantities for a weak 
and a strong first order phase transition. In some limiting cases their behavior can be 
studied analytically. For the discussion we will frequently consider the flow equations for 
the quartic 'couplings' Ai(p), \2{p) eqs. (6.8), (6.9) and for the minimum k eq. (6.15). 

In flg. 6, 7 we follow the flow of the dimensionless renormalized minimum k and 
the radial mass term fh? = 2kXi{k) in comparison to their dimensionful counterparts 
PoR — kn and mji = k'^m^ in units of the momentum scale A. We also consider the 
dimensionless renormalized mass term 7ri| = kA2(k) corresponding to the curvature of the 
potential in the direction of the second invariant f. The height of the potential barrier 
Usik) = k^Uk{pB) with u'f.{pB) = 0, < Pb < k, and the height of the outer minimum 
Uo{k) = k^Uk{n) is also displayed and will be discussed in section 9. Fig. 6 shows these 
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t = ln(A;/A) 

Figure 6: Scale dependence of the dimensionless renormalized masses m^, ml, minimum k and dimen- 

sionful counterparts = k'^fh^, p^R = kn in units of A. Wc also show [/^(fc) and Uaik), the value of 
the potential at the top of the potential barrier and at the minimum poij, respectively. The short distance 
parameters are Aia = 2, A2A = 0.1 and Ska = 0. 

quantities as a function of t — ln(/c/A) for Aia = 2, A2A = 0.1. One observes that the 
flow can be separated into two parts. The first part ranging from t = to t ~ —6 is 
characterized by k ~ const and small 7712- It is instructive to consider what happens in 
the case ml = 0. In this case A2 = and the flow is governed by the Wilson-Fisher flxed 
point of the 0(8) symmetric theory. At the corresponding second order phase transition 
the evolution of Uk leads to the scahng solution of (5.15) which obtains for duk/dt — 0. As 
a consequence becomes a fc-independent function that takes on constant (fixed point) 
values [27, 22]. In particular, the minimum k of the potential takes on its fixed point 
value K,{k) = n^. The fixed point is not attractive in the U{2) x U (2) symmetric theory 
and A2A is an additional relevant parameter for the system. For small A2 the evolution is 
governed by an anomalous dimension d\2/dt — AX^ with ^4 < 0, leading to the increasing 
rfil as k is lowered. 

The system exhibits scaling behavior only for sufficiently small A2. As m| increases 
the quartic coupling Ai and therefore the radial mass term iv? is driven to smaller values 
as can be observed from fig. 6. For nonvanishing A2 the corresponding qualitative change 
in the flow equation (6.8) for Ai is the occurance of a term ~ A|. It allows to drive Ai 
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to negative values in a certain range oi p < k and, therefore, to create a potential barrier 
inducing a first order phase transition. We observe from the plot that at — 9.5 a second 
minimum arises {Ub ^ 0)- The corresponding value of A; = Ae* = ^2 sets a characteristic 
scale for the first order phase transition. Below this scale the dimensionless, renormalized 
quantities approximately scale according to their canonical dimension. The dimensionful 
quantities like Pqr or mj^ show only a weak scale dependence in this range. In contrast to 
the above example of a weak first order phase transition with characteristic renormalized 
masses much smaller than A fig. 7 shows the fiow of the corresponding quantities for a 
strong first order phase transition. The short distance parameters employed are Aia = 0.1, 
A2A = 2. Here the range with k ~ const is almost absent and one observes no approximate 
scaling behavior. 




-2 -1.5 -1 -0.5 



t = ln{k/A) 
Figure 7: Same as fig. 6 for Aia =0.1 and A2A = 2. 

With the above examples it becomes easy to understand the phase structure presented 
in section 4. For the given curves of figs. 4 and 5 we distinguished between the range 
A2A/A1A <S 1 and A2A/A1A ^ 1 to denote the weak and the strong first order region. For 
A2A/A1A ^ 1 the initial renormalization group flow is dominated by the Wilson-Fisher 
fixed point of the 0(8) symmetric theory. In this range the irrelevant couplings are driven 
close to the fixed point for some 'time' \t\ = — ln(A;/A), loosing their memory on the 
initial conditions given by the short distance potential u\. As a consequence we are able 
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to observe universal behavior as is demonstrated in fig. 5. 

To discuss the case A2A/A1A ^ 1 we consider the flow equations for the couphngs at 
the minimum k 7^ of the potential given by (6.15) and (6.16) with (6.8), (6.9). In the 
limit of an infinite mass term = k\2{k) 00 the /9-functions for Ai(k) and k, become 
independent from A2(/t) due to the 'threshold' functions, with 1^(k\2) ~ {k\2)^''"''^^^ for 
large k,X2{k). As a consequence and equal the /3-functions for an 0(5) symmetric 
model. We argue in the following that in this large coupling limit fluctuations of massless 
Goldstone bosons lead to an attractive flxed point for X^ii^). We take the flow equation 
(6.16), (6.9) for A2(k) keeping only terms with positive canonical mass dimension for a 
qualitative discussion. (This amounts to the approximation X^I^\k) = \2^\k) = for 
n > 1.) To be explicit, one may consider the case for given Aia = 2. The critical cutoff 
value for the potential minimum is Ka — 0.2 for A2A ^ 1- For kA2(k) 3> 1 and taking 
?7 ~ the /3-function for X2{k) is to a good approximation given by {d — 3) 

dA2(«) , , X , ^_ , ,.2,3 



dt 



^-X2{K) + 2v,{X2{K)yil(0). (7.1) 



The second term on the r.h.s. of eq. (7.1) is due to massless Goldstone modes which 
give the dominant contribution in the considered range. The solution of (7.1) implies an 
attractive flxed point for A2(fi;) with a value 

Starting from A2A one flnds for the 'time' \t\ necessary to reach a given X2{k) > A2*(k) 



A.(.) (1 - ^ 



This converges to a flnite value for A2A — ^ 00. The further evolution therefore becomes 
insensitive to the initial value for A2A in the large coupling limit. The flow of Xi{k) and 
K is not affected by the initial running of A2(k) and quantities like Apofi/A or mR/ApoR 
become independent of A2A ff the coupling is sufficiently large. This qualitative discussion 
is confirmed by the numerical solution of the full set of equations presented in figs. 4 and 
5. For the fixed point value we obtain A2*(«;) = 38.02. We point out that an analogous 
discussion for the large coupling region of Aia cannot be made. This can be seen by 
considering the mass term at the origin of the short distance potential (3.6) given by 
m'^(0,0) = — kaAia- Due to the pole of l^{uj,rj) at a; = — 1 for n > 1/2 (cf. appendix A) 
one obtains the constraint 

«aAia < 1 . (7.4) 

In the limit Aia 00 the mass term 2kaAia at the minimum k of the potential at the 
critical temperature therefore remains finite. 
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8 Scaling equation of state for weak first order phase 
transitions 



We presented in section 4 some characteristic quantities for the effective average potential 
which become universal at the phase transition for a sufficiently small quartic coupling 
A2A = A2A/A of the short distance potential U\ (3.6). The aim of this section is to 
generalize this observation and to find a universal scaling form of the equation of state for 
weak first order phase transitions. The equation of state relates the derivative of the free 
energy U — limfc^o t^fc to an external source, dU/dip — j. Here the derivative has to be 
evaluated in the outer convex region of the potential. For instance, for the meson model of 
strong interactions the source j is proportional to the average quark mass [8, 36] and the 
equation of state permits to study the quark mass dependence of properties of the chiral 
phase transition. We will compute the equation of state for a nonzero coarse graining 
scale k. It therefore contains information for quantities like the 'classical' bubble surface 
tension in the context of Langer's theory of bubble formation which will be discussed in 
section 9. 

In three dimensions the U{2) x U (2) symmetric model exhibits a second order phase 
transition in the limit of a vanishing quartic coupling A2A due to an enhanced 0(8) 
symmetry. In this case there is no scale present in the theory at the critical temperature. 
In the vicinity of the critical temperature (small \Sk,\\ ~ |Tc — r|) and for small enough 
A2A one therefore expects a scaling behavior of the effective average potential Uk and 
accordingly a universal scaling form of the equation of state. At the second order phase 
transition in the 0(8) symmetric model there are only two independent scales that can 
be related to the deviation from the critical temperature and to the external source or 
(p. As a consequence the properly rescaled potential C^/p|. or [//p^^+^^Z^ (with the usual 
critical exponent 5) can only depend on one dimensionless ratio. A possible set of variables 
to represent the two independent scales are the renormalized minimum of the potential 
^OR = (poij/2)^^^ (or the renormalized mass for the symmetric phase) and the renormalized 
field cpn — (pr/2)^/^. The rescaled potential will then only depend on the scahng variable 
z = i^r/^or [15]. Another possible choice is the Widom scaling variable x — —SK\/ip^^^ 
[38]. In the U{2) x U{2) symmetric theory A2A is an additional relevant parameter which 
renders the phase transition first order and introduces a new scale, e.g. the nonvanishing 
value for the jump in the renormalized order parameter Alpqr = (Apor/2)^''^ at the critical 
temperature or 5ka — 0. In the universal range we therefore observe three independent 
scales and the scaling form of the equation of state will depend on two dimensionless 
ratios. The rescaled potential U/ip^j^ can then be written as a universal function G 



U 



G{z,v) 




which depends on the two scaling variables 



^R 



(8.2) 



z = 



V — 



^QR 



<Pqr 
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The relation (8.1) is the scahng form of the equation of state we are looking for. At a 
second order phase transition the variable v vanishes identically and G{z,0) describes the 
scaling equation of state for the model with 0(8) symmetry [15]. The variable v accounts 
for the additional scale present at the first order phase transition. We note that z = 1 
corresponds to a vanishing source and G{l,v) describes the temperature dependence of 
the free energy for j = 0. In this case v = 1 denotes the critical temperature whereas 
for T < Tc one has v < 1. Accordingly v > 1 obtains for T > Tc and (poR describes here 
the local minimum corresponding to the metastable ordered phase. The function G{z, 1) 
accounts for the dependence of the free energy on j for T = T^. 

We consider the scaling form (8.1) of the equation of state for a nonzero coarse graining 

1/2 

scale k. The renormalized field is given by (pn = ip. We pointed out in section 4 that 
there is a characteristic scale A;2 for the first order phase transition where the second local 
minimum of the effective average potential appears. For weak first order phase transitions 
one finds poR ~ ^2- To observe the scaling form of the equation of state the infrared cutoff 
k has to run below k2 with k -C /c2- For the scale kf defined in eq. (4.11) we observe 
universal behavior to high accuracy (cf. fig. 5 and the corresponding universal ratios in 
table 1 for small A2a/Aia). The result for the universal function UkJ^QR = Gkf{z,v) is 
presented in fig. 8. For v — 1 one has (poR{kf) — Aiponikf) which denotes the critical 
temperature. Accordingly v > 1 denotes temperatures above and v < 1 temperatures 
below the critical temperature. One observes that Gkf{z, 1) shows two almost degenerate 
minima. (They become exactly degenerate in the limit /c — > 0). For the given examples 
V — 1.18, 1.07 the minimum at the origin becomes the absolute minimum and the system 
is in the symmetric phase. In contrast, for v — 0.90, 0.74 the absolute minimum is located 
at ^ = 1 which characterizes the spontaneously broken phase. For small enough v the 
local minimum at the origin vanishes. 

We explicitly verified that the universal function Gkj depends only on the scaling 
variables z and v by choosing various values for 6k,a and for the quartic couplings of 
the short distance potential, Aia and A2A. In section 4 we observed that the model 
shows universal behavior for a certain range of the parameter space. For given Aia and 
small enough A2A one always observes universal behavior. For Aia = 0.1, 2 and 4 it is 
demonstrated that (approximate) universality holds for A2a/Aia^ 1/2 (cf. fig. 5). For Aia 
around 2 one observes from figs. 4, 5 and table 1 that the system is to a good accuracy 
described by its universal properties for even larger values of A2a- The corresponding 
phase transitions cannot be considered as particularly weak first order. The universal 
function Gkf therefore accounts for a quite large range of the parameter space. 

We emphasize that the universal form of the effective potential given in fig. 8 depends 
on the scale kf where the integration of the fiow equations is stopped (cf. eq. (4.11)). 
A different prescription for kf will, in general, lead to a different form of the effective 
potential. We may interpret this as a scheme dependence describing the effect of different 
coarse graining procedures. This is fundamentally different from nonuniversal corrections 
since Gk^ is completely independent of details of the short distance or classical action 
and in this sense universal. A more quantitative discussion of this scheme dependence 
will be presented in the next section. We note that fluctuations on scales k < kf do not 
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Figure 8: Universal shape of the coarse grained potential {k ^ kf) as a function of the scaling variable 
z = (pn/ipm = {pr/porY''^ for different values oi v = Afon/foR = {^por/ PorV^^ ■ The employed 
values for v are (1) v — 1.18, (2) v = 1.07, (3) w = 1, (4) w = 0.90, (5) v — 0.74. For vanishing sources 
one has z = I. In this case v = 1 denotes the critical temperature Tc- Similarly v > 1 corresponds to 
T > Tc with ifoji denoting the minimum in the metastable ordered phase. 

influence substantially the location of the minima of the coarse grained potential and the 
form of Uk{(pR) for Lpf^ > ipQn remains almost /c-independent. Here dU^Jdip — j{kf) with 
j{kf ) ~ limfc^oj(^) = j-^ 

We have established the scaling equation of state using the renormalized variables ipR, 
(foR and A(foR. Alternatively the scaling equation of state can be presented in terms of 
the short distance parameters 5k\, A2A and the unrenormalized fleld (p. For the 0{8) 
symmetric model (A2A = 0) this is known as the Widom scaling form of the equation of 
state [38] and in this case the relation between (pR, Pqr and (f, 6k,\ is solely determined 
by critical exponents and amplitudes [15]. For the U{2) x U{2) symmetric model the 
dependence of (fR, ipoR and A<^oi? on the parameters Sk\, A2A and the unrenormalized 
field (fi can be expressed in terms of scaling functions. In general, only for limiting cases as 
5k\ = or A2A = the relation is determined by critical exponents and amplitudes. We 
will consider here these limits. The computation of the scaling functions for the general 
case as well as the corresponding (generalized) Widom scaling form of the equation of 
state is the subject of a separate work [39]. 

^The role of massless Goldstone boson fluctuations for the universal form of the effective average 
potential in the limit A; — > has been discussed previously for the 0(8) symmetric model [15]. 
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We consider the renormalized minimum (^qr in two limits which are denoted by 
^VoR = ^or{Sk\ — 0) and ip^ji = (Por{X2A — 0). The behavior of A(/7or is described 
in terms of the exponent 9 according to eq. (4.7), 

Av^oi? - (A2a)'/', = 1.93. (8.3) 

The dependence of the minimum ip^^^ of the 0(8) symmetric potential on the temperature 
is characterized by the critical exponent u, 

foR ~ {SkaT^', y = 0.882. (8.4) 

The exponent v for the 0(8) symmetric model is determined analogously to Q as de- 
scribed in section 4 ^. We can also introduce a critical exponent ^ for the jump of the 
unrenormalized order parameter 

A(^o-(A2a)^, C = 0.988. (8.5) 

With 

~ (5«a)^, /9 = 0.451 (8.6) 
it is related to Q and v by the additional index relation 

^^^ = 05. (8.7) 

We have verified this numerically. For the case 5ka = A2A = one obtains 

J ~ (8.8) 

The exponent 6 is related to the anomalous dimension 77 via the usual index relation 
5 = (5 — 7;)/(l + T]). From the scaling solution of eq. (6.21) we obtain r] = 0.0224. 
With the help of the above relations one immediately verifies that for A2A = 

zr.i-x)-'^, v^O (8.9) 

and for 5k\ — 

z^y-'^, v^l. (8.10) 
Here we have used the Widom scaling variable x and the new scaling variable y given by 



'^For the 0{8) symmetric model (A2A = 0) we consider the minimum tpg^j at A: = 0. 
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9 Coarse graining and Langer's theory of bubble for- 
mation 



The coarse grained effective potential C4 results from the integration of fluctuations with 
momenta larger than ~ k. It is a nonconvex function whereas the standard effective 
potential U = limfc^o Uk has to be convex by its definition as a Legendre transform. 
The difference between U and Uk is due to fluctuations with characteristic length scales 
larger than the inverse coarse graining scale ~ The role of these fluctuations for 

the approach to convexity has been estabhshed explicitly [26]. The study of flrst order 
phase transitions usually relies on the nonconvex 'part' of the potential. As an illustration 
we consider the change of the system from the high temperature to the low temperature 
phase by bubble nucleation as described by Langer's theory [17]. On the one hand, the 
approach rehes on the deflnition of a suitable coarse grained free energy with a coarse 
graining scale k and, on the other hand, a saddle point approximation for the treatment of 
fluctTiations around the 'critical bubble' is employed. The problem is therefore separated 
into two parts: One part concerns the treatment of fluctuations with momenta ^ k'^ 
which are included in the coarse grained free energy. The second part deals with an 
estimate of fluctuations around the bubble for which only fluctuations with momenta 
smaller than k must be considered. These issues will be discussed in this section in 
a quantitative way. In particular, we will give a criterion for the validity of Langer's 
formula. 

One may consider a system that starts at some high temperature T > and in- 
vestigate what happens as T is lowered as a function of time as for example during the 
evolution of the early universe. For large enough temperature the origin of the potential 
(</? = 0) is the only minimum and the system is therefore originally in the symmetric 
phase. As T approaches a second local minimum develops at (/9o > 0. This becomes 
the absolute minimum below Tc. Nevertheless, the potential barrier prevents the system 
to change smoothly to the ordered phase. For a short while where T is in the vicinity 
of Tc but below Tc the system remains therefore in a state with higher energy density as 
compared to the state corresponding to the absolute minimum away from the origin. This 
is the so-called 'false vacuum' in high energy physics or the metastable state in statistical 
physics. Such a state is unstable with respect to fluctuations which penetrate or cross 
the barrier. The picture is familiar from the condensation of vapor. The false vacuum 
corresponds to the supercooled vapor phase and the true vacuum to the fluid phase. Bub- 
bles of the true vacuum (droplets) occur through thermal or quantum fluctuations^. If a 
bubble is large enough so that the decrease in volume energy exceeds the surface energy 
it will grow. The phase transition is completed once the whole space is filled with the 
true vacuum. On the other hand, small bubbles shrink due to the surface tension. The 
critical bubble is just large enough that it does not shrink. To be explicit we consider 
a spherical bubble where the bubble wall with 'thickness' A is thin as compared to the 
bubble radius i?, i.e. A ^ i?. In leading order the coarse grained free energy Fj. for such 

^In the real world the condensation of vapor is triggered by impurities but this is not the issue here. 
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a bubble configuration can be decomposed in a volume and a surface term [40, 41] 

rf^-^Rh + ATTRW (9.1) 

o 

In the thin wall approximation one obtains for the surface tension ak in our conventions 

<po 

ak^2 1 dipy/2ZkUk{<p). (9.2) 



For the difference in the free energy density e one has 

e = U{0) - U{ipo) = - lim Uo{k). (9.3) 

fe— >o 

We include in e fluctuations with arbitrarily small momenta. In contrast, the long wave- 
length contributions to 0"^ are effectively cut off by the characteristic length scale of the 
bubble surface and are described by the 'fluctuation determinant' (cf. eq. (9.4)). The 
determination of a suitable coarse graining scale k for the computation of ak is discussed 
below. 

The critical bubble maximizes T^^^ with respect to the radius. It minimizes the coarse 
grained free energy with respect to other deformations because the spherical form is 
energetically favorable. The critical bubble therefore represents a saddle point in the space 
of possible 'bubble configurations'. In hanger's theory of bubble formation one considers 
a saddle point expansion around the critical bubble. There is exactly one negative mode 
that corresponds to the shrinking or growth of the bubble and there are infinitely many 
positive modes (there are also translational zero modes) ^. The bubble nucleation rate T, 
which describes the probability per unit volume per unit time for the transition to the 
new vacuum, can be written in the form [40, 41] 

f = A exp(-rl°)[^r]) = A, exp {-'-ffj (9.4) 

where F^^"* is evaluated for ipf'^ corresponding to the critical bubble and approximated 
by (9.1). The exponential term with the coarse grained free energy f[.°^ denotes the 
lowest order or classical contribution. The prefactor contains several factors that 
depend on the details of the system under investigation. In particular, A^ accounts for 
the contribution to the free energy from the fluctuations with momenta smaller than 
k. It depends on k through the effective ultraviolet cutoff for these fluctuations which is 
present since fluctuations with momenta larger than k are already included in F^°'*[(^[°'*] 

^Langer's theory is not restricted to the thin wall approximation which is considered here for simplicity. 
In particular, the property of the critical bubble to represent a saddle point with exactly one negative 
mode is independent from the thin wall approximation. 

^°The effective average action [18] also provides the formal tool how the ultraviolet cutoff ~ fc is 
implemented in the remaining functional integral for large length scale fluctuations. 
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Langer's formula for bubble nucleation amounts essentially to a perturbative one loop 
estimate of ^4^. 

For a determination of a useful choice of k it is convenient to place the discussion in 

a more general context which does not rely on the thin wall approximation or a saddle 
point approximation. What one is finally interested in is the free energy T[lpi,] for bubble 
configurations of a given shape. The 'true critical bubble' ip^ corresponds to a saddle point 
in the space of 'bubble configurations' which are characterized by boundary conditions 
connecting the false and the true vacuum. The nucleation rate is then proportional 
exp— r[(pg]. The coarse graining can be seen as a convenient strategy to evaluate r[(p;,] 
by separating contributions from different momentum scales. We propose to choose the 
coarse graining scale k somewhat above but in the vicinity of the inverse thickness 
of the bubble wall. We will argue below that in this case the corrections to the effective 
surface tension from the prefactor Ak should be best accessible. 

In fact, we can write T = B exp — r[(^^] where B contains dynamical factors and T[ipl] 
does not include contributions from fluctuations of the negative mode and the translational 
modes present for the critical bubble. The prefactor in eq. (9.4) can then be written as 

Ak ^ B exp -(5k + Vk) (9.5) 

where 

6k = TK] - r^K] , vk = r^K] - — ^ . (9.6) 

The term f]k includes the difference between the true critical bubble and the configuration 
used to estimate ak as well as a correction term to e to be discussed below. We first 
concentrate on 5k which describes the difference between the free energy and the coarse 
grained free energy for the critical bubble. As mentioned above this is due to fluctuations 
with momenta q'^^k'^ and incorporates the dominant ^-dependence of A^. Since the 
bubble provides for inherent effective infrared cutoff scales ~ R^^ or the contribution 
6k is both infrared and ultraviolet finite. The larger k, the more fiuctuations are included 
in 5k and from this point of view one wants to take k as low as possible. On the other hand, 
k should not be taken smaller than if the approximation used for a computation of 
Tk relies on almost constant field configurations rather than real bubbles, as is usually the 
case. Only for k sufficiently large compared to A^^ the difference between an evaluation 
of the potential and kinetic terms in for almost constant field configurations (e.g. by 
a derivative expansion) rather than for bubbles remains small. In this way the technique 
of course graining combines a relatively simple treatment of the modes with ^ k'^ for 
which the detailed properties of the bubble are irrelevant with an estimate of fiuctuations 
around the bubble for which the short distance physics (g^ ^ k"^) needs not to be considered 
anymore. It is clear that k is only a technical construct and for physical quantities 
the fc-dependence of 5k and must cancel. More precisely, this concerns the sum 
5k + Vk + 167r(j|/3e^ . For thin wall bubbles the most important contribution to r]k is 
easily identified: By our definition of e we have included contributions from fiuctuations 
with length scales ^R. They should not be present in the effective action for a bubble 
with finite radius. Therefore rjk contains a correction term (167r/3)cr|(e~^(i?) — e~^) which 
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replaces effectively e by e(i?) in eq. (9.4). We can evaluate e{R) in terms of the coarse 
grained free energy at a scale kji 



e{R) C/,,(0) - Uk.i^o) 



Uo{kR) , 




(9.7) 



For A <^ R one should not confound kn with the coarse graining scale k since one has 
the inequality 



Only for A ~ i? the clear separation between kji and k disappears. We note that at the 
critical temperature one has i? — > cxo and therefore e(/?) = e . 

Since 0"^ enters the nucleation rate (9.4) exponentially even small changes with k will 
have large effects. If one finds a strong dependence of ak on the coarse graining scale k 
this is only compatible with a large contribution from the higher orders of the saddle point 
expansion. The fc-dcpendence of <Tk therefore gives direct information about the validity 
of hanger's formula. We find a strong scale dependence of if the phase transition is 
characterized by large effective dimensionless couphngs ^{k). A weak scale dependence is 
observed for small effective couplings. This gives a very consistent picture: The validity 
of the saddle point approximation typically requires small dimensionless couplings. In 
this case also the details of the coarse graining are not of crucial importance within an 
appropriate range of k. The remaining part of this section provides a detailed quantitative 
discussion. 

We consider the dependence of the effective average potential 17^ and the surface 
tension ak on the coarse graining scale k near the critical temperature T^. for three examples 
in detail. They arc distingTiishcd by different choices for the quartic couplings Aia and 
A2A of the short distance potential Ua given by eq. (3.6). The choice Aia/A = 0.1 and 
A2A/A = 2 corresponds to a strong first order phase transition with renormalized masses 
not much smaller than the cutoff scale A. The renormahzed couphngs will turn out small 
enough such that the notion of a coarse grained potential Uk and a surface tension can 
be used without detailed information on the coarse graining scale within a certain range 
of k. In contrast we give two examples where the dependence of Uk and (Tk on the coarse 
graining scale becomes of crucial importance. The choice Aia/ A = 2 and A2A/A = 0.1 
leads to a weak first order phase transition with small renormalized masses and the system 
shows universal behavior (cf. sect. 4). For Aia/A = 4 and A2A/A = 70 one observes a 
relatively strong first order phase transition. Nevertheless for both examples the coarse 
grained potential and the surface tension show a similarly high sensitivity on the scale k. 

In figs. 9, 10 the scale dependence oiUk is shown for a fixed temperature in the vicinity 
of T = Tc or 5ka = 0. We plot Uk in units of A^ as a function of ^Pr/A^^'^ = (pr/2A)1/2 ^ 
(Zfcp/2A)^/^ for different values of k. Each curve differs in A; = Ae* by At — 1/18 and the 
first curve to the left with a negative curvature at the origin corresponds to t ~ —0.40 
for fig. 9 (t ~ -9.3 for fig. 10). For > t;^ - 0.40 (0 > t;^ - 9.3) there is only one 
minimum of the potential away from the origin which lies below the plotted range in 
fig. 9 (10). Lowering k results in a successive inclusion of fiuctuations on larger length 



kR < -r ^ . 
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Figure 9: The coarse grained potential in dependence on the coarse graining scale k for fixed (almost 
critical) temperature. The example corresponds to Aia/A = 0.1, A2a/A = 2. 

scales and one observes the appearance of a second local minimum at the origin of the 
potential. The barrier between both minima increases and reaches its maximum value. 
As k is lowered further the potential barrier starts to decrease whereas the location of the 
minima becomes almost /c-independent and degenerate in height. In this region also the 
outer convex part of the potential shows an almost stable profile due to the decoupling 
of the massive modes. We stop the integration at the scale k — kf given by eq. (4.11). 
The coarse grained effective potential Uk with k = kf essentially includes all fluctuations 
with squared momenta larger than the scale given by the curvature at the top of 

the potential barrier (cf. eq. (4.10)). Successive inclusion of fluctuations with momenta 
smaller than kf would result in a further decrease of the potential barrier. The flattening 
of the barrier is induced by the pole structure of the 'threshold' functions l^i^) appearing 
in the flow equations for the couplings (e.g. iKou) appearing in eq. (6.8) exhibits a pole at 
u = —1 according to /f (c^) ~ {to + 1)~^/^ for to near —1 (cf. appendix A)). The argument 
uj corresponds to dimensionless mass terms as U'^^^pRj/k'^ or {U'j^{pR) + 2prU'^{pr)) / k"^ . In 
the nonconvex region of the potential the curvature is negative. Since the pole cannot be 
crossed, negative Ul or U'^ must go to zero with A;^ and as a consequence the potential 
barrier flattens. For the scalar 'vector' model the approach to convexity in the limit /c — 
has been studied analytically previously [26]. Here we are interested in the potential for 
a nonzero coarse graining scale k. 
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Figure 10: The coarse grained potential for Aia/A = 2, A2A/A = 0.1. 



The most significant difference between fig. 9 and fig. 10 is the /^-dependence of the 
potential barrier in a region where the minima become degenerate and almost independent 
of k. Fig. 9 shows a barrier with a weak scale dependence in the range of k where the 
location of the minima stabilizes. In contrast, in fig. 10 one observes a barrier with a 
strong scale dependence in this region. Figs. 6 and 7 (cf. sect. 7) show the relative height 
Uo{k) = f/fc((/9o) — Uk{0) between the two local minima and the height of the potential 
barrier Usik) = Uk{^B) — Uk{0) {{dUk/ d(p){(pB) = 0, < < (po) as a function of 
t — ln(A;/A). Accordingly one observes that for Aia/A — 0.1, A2A/A = 2 the top of the 
potential barrier shows a weak scale dependence in a region of k where UQ{k) is small 
whereas for Aia/A = 2, A2A/A = 0.1 the top of the potential barrier depends strongly 
on the coarse graining scale in this region. The surface tension is displayed in fig. 
11 and shows a corresponding behavior. Here we consider ak also for the short distance 
parameters Aia/A = 4, A2A/A = 70. In fig. 11 the surface tension Ok is normalized to 
(Jmax (w/A' = 1-67 X 10-2(8.41 X 10-ii)(1.01 X 10-=^) for Aia/A = 0.1(2)(4), Asa/A = 
2(0.1)(70)) and given as a function of \n{k/kf)}^ For Aia/A = 0.1, A2A/A = 2 the curve 
exhibits a small curvature around its maximum and fTmax — <^k=kf One observes for 
the second and the third example a comparably large curvature around cxmax and ak=kf 
becomes considerably smaller than the maximum. One may consider what happens if 



^^The integration according to eq. (9.2) is performed between the two zeros = and ip = i^'q ^(fio of 
Uk. 
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\n{k/kf) 

Figure 11: The normalized surface tension (Tfe/(Tmax as a function of ln(fc/fc/). The short distance 
parameters are (1) Aia/A = 0.1, A2a/A = 2, (2) Aia/A = 2, W/A = 0.1, (3) Aia/A = 4, W/A = 70. 



this scale dependence is not taken into account correctly. If one takes the difference in at 
from its maximum value to (Tk=kf as a rough measure for its uncertainty this is about 30% 
for the last example. Entering exponentially in eq. (9.4) this would lead to tremendous 
errors. 

In order to quantify the differences between the three examples we have displayed 
some characteristic quantities in table 2. The renormalized couplings 

AlR=[/fe,(pOR), X2R^^drUk,{poR) (9.9) 

are normalized to the mass term 



m 



R 



(2porAik)V^ 



In addition we give the mass term 



nT'2R = {PoR>^2Ry^^ 



(9.10) 



(9.11) 



corresponding to the curvature of the potential in the direction of the second invariant 
T. In comparison with figs. 9-11 one observes that the smaller the effective couplings 
the weaker the scale dependence of Uk and (7^. In particular, a reasonably weak scale 
dependence of Uk and ak requires 

XiR 1 ml 



m' 



R 



'OR 



< 1 . 



(9.12) 
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^2 A 




A2_R 










T 


~T 


m% 








A 


A 


0.1 


2 


0.228 


8.26 


0.235 


1.55 X 10"^ 


6.62 X 10-1 


1.011 X 10-1 


2 


0.1 


0.845 


15.0 


0.335 


2.04 X 10-5 


6.10 X 10-5 


1.145 X 10-5 


4 


70 


0.980 


16.8 


0.341 


6.96 X 10-2 


2.04 X 10-1 


3.781 X 10-2 



Table 2: The effective dimensionless couplings Xm/m'^ and A2i{/m^. The couphngs and the mass 
terms m^, are evaluated at the scale kf and 5k\ = 0. 

This establishes a quantitative criterion for the range where Langer's theory can be used 
without paying too much attention to the precise definition of the coarse graining. Com- 
parison with fig. 5 (cf. sect. 4) shows that this condition is not realized for the range 
of couplings leading to universal behavior and for large Aia/A. The only region where 
the saddle point expansion is expected to converge reasonably well is for small Aia/A 
and large A2a/A. The second and the third example given in fig. 11, which exhibit a 
strong /c-dependence, show similar values for the effective couplings. More precisely, for 
the relatively strong phase transition with slightly larger effective couplings one observes 
an increased scale dependence as compared to the weak phase transition. In table 2 also 
the renormalized masses in units of A which indicate the strength of the phase transition 
and kf/A are presented. 

In summary, we have shown that the coarse grained free energy cannot be defined 
without detailed information on the coarse graining scale k unless the effective dimen- 
sionless couplings are small. Only for small couplings we observe a weak /c-dcpcndcnce of 
the surface tension in a range where the location of the minima of the potential remains 
almost fixed. There is a close relation between the dependence of the coarse grained free 
energy on the coarse graining scale and the reliability of the saddle point approximation in 
Langer's theory of bubble nucleation. For a strong fc-dependence of Uk a small variation in 
the coarse graining scale can induce large changes in the predicted nucleation rate in low- 
est order in a saddle point approximation. In this case the /c-dependence of the prefactor 
Ak has also to be computed. For strong dimensionless couphngs a realistic estimate of the 
nucleation rate therefore needs the capability to compute In Ak with the same accuracy as 
167ia'l/3e^ and a check of the cancelation of the fc-dependence in the combined expression 
(9.4). Our observation that the details of the coarse graining prescription become less 
important in the case of small dimensionless couplings is consistent with the fact that 
typically small couplings are needed for a reliable saddle point approximation for Ak. For 
the electroweak high temperature phase transition a small /c-dependence of ak is found 
for a small mass Mh of the Higgs scalar whereas for Mh near the W-boson mass the 
picture resembles our fig. 10 [42]. This corresponds to the observation [43, 44, 45] that 
the saddle point approximation around the critical bubble converges well only for a small 
enough mass of the Higgs scalar. 
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10 Conclusions 



We have presented in this paper a detailed investigation of the phase transition in three 
dimensional models for complex 2x2 matrices. They are characterized by two quartic 
couplings AiA and A2a- In the limit Aia — oo, X2A 00 this also covers the model of 
unitary matrices. The picture arising from this study is unambiguous: 

(1) One observes two symmetry breaking patterns for A2A > and A2A < respectively. 
The case A2A = denotes the boundary between the two phases with different symmetry 
breaking patterns. In this special case the theory exhibits an enhanced 0(8) symmetry. 
The phase transition is always first order for the investigated symmetry breaking U (2) x 
U{2) U{2) (A2A > 0). For A2A ~ the 0(8) symmetric Heisenberg model is recovered 
and one finds a second order phase transition. 

(2) The strength of the phase transition depends on the size of the classical quartic 
couplings Aia/A and A2A/A. They describe the short distance or classical action at a 
momentum scale A. The strength of the transition can be parametrized by m^/A with 
rrfji a characteristic inverse correlation length at the critical temperature. For fixed A2A 
the strength of the transition decreases with increasing Aia- This is analogous to the 
Coleman- Weinberg effect in four dimensions. 

(3) For a wide range of classical couplings the critical behavior near the phase transition 
is universal. This means that it becomes largely independent of the details of the classical 
action once everything is expressed in terms of the relevant renormalized parameters. In 
particular, characteristic ratios like m^/ApoR (critical inverse correlation length in the 
ordered phase over discontinuity in the order parameter) or m§^/Apo_R (same for the 
disordered phase) arc not influenced by the addition of new terms in the classical action 
as far as the symmetries are respected. 

(4) The range of short distance parameters Aia, A2A for which the phase transition 
exhibits universal behavior is not only determined by the strength of the phase transition 
as measured by rrf^/ K. For a given Aia/A and small enough A2A/A one always observes 
universal behavior. In the range of small Aia/A the essential criterion for universal be- 
havior is given by the size of A2A/A1A, with approximate universality for A2A < Aia- For 
strong couplings universality extends to larger A2A/A1A and occurs for much larger m^/A 
(cf. table 1). 

(5) We have investigated how various characteristic quantities like the discontinuity 
in the order parameter Apo or the corresponding renormalized quantity Apoi? or criti- 
cal correlation lengths depend on the classical parameters. In particular, at the critical 
temperature one finds universal critical exponents for not too large A2A, 

Apoi? ~ (W, ^^ = 1.93, 

Apo ~ (A2a)'^, C = 0.988 . 

These exponents are related by a scaling relation to the critical correlation length and 
order parameter exponents v and j3 of the 0(8) symmetric Heisenberg model according 
to e/C =_T^/I3 = 1.95 (z/ = 0.882, fS = 0.451 in our calculation for Asa = 0). Small 
values of A2A can be associated with a perturbation of the 0(8) symmetric model and 9, ( 
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are related to the corresponding crossover exponents. On the other hand, ApoR (Apo) 
becomes independent of A2A in the infinite couphng hmit. 

(6) We have computed the universal equation of state. The equation of state relates 
the derivative of the free energy U to an external source, dU/dip = j. From there 
one can extract universal ratios e.g. for the jump in the order parameter [ApQ^/nf^ = 
0.592) or for the ratios of critical correlation lengths in the disordered (symmetric) and 
ordered (spontaneously broken) phase {rriQ^/m'j^ — 0.746). It specifies critical couplings 
(Aiij/m^ = 0.845, XiR/rrfji — 15.0). The universal behavior of the potential for large field 
arguments p/j ^ po_R is [/ ~ pfj ~ p3/(i+»7) provided p/j is sufficiently small as compared 
to A. Here the critical exponent t] which characterizes the dependence of the potential on 
the unrenormalized field p is found to be 77 = 0.022. For large p the universal equation 
of state equals the one for the 0(8) symmetric Heisenberg model and r] specifies the 
anomalous dimension or the critical exponent 5 — {b — r])/{l-\-r]). The equation of state 
is computed for a nonzero coarse graining scale k. It therefore contains information for 
quantities like the 'classical' bubble surface tension in the context of Langer's theory of 
bubble formation. 

(7) We have investigated the dependence of the coarse grained effective potential 
Uk{p) and the 'classical' surface tension Uk on the coarse graining scale k with special 
emphasis on the question of the validity of Langer's theory of bubble formation. Wc 
find a strong scale dependence of f/^ and if the phase transition is characterized by 
large dimensionless couplings. A weak scale dependence is observed for small effective 
couphngs. There is a close relation between the dependence of the coarse grained free 
energy on the coarse graining scale and the reliability of the saddle point approximation 
in Langer's theory of bubble nucleation. A strong /c-dependence of ak is only compatible 
with a large contribution from the higher orders of the expansion. We obtain a very 
consistent picture: The validity of the saddle point approximation typically requires small 
dimensionless couplings. In this case also the details of the coarse graining are not of 
crucial importance within an appropriate range of k. The quantitative criterion for the 
validity of Langer's formula is in o\ii case Ai/,>/m^ ^ 1. 

(8) Our method is not restricted to the study of the universal behavior. We can 
compute the effective potential for arbitrary values of the initial parameters and have 
done this for particular examples. 

The uncertainties of our results induced by the numerical integration of the fiow equa- 
tions are well under control and small. They are negligible compared to the expected 
error induced by our truncation. For a significant improvement of our treatment one 
would have to include higher order terms in the derivative expansion employed for the 
effective average action. For weak first order or second order phase transitions we expect 
the error to be related to the anomalous dimension 77 = 0.022. For the special case of the 
enhanced 0(8) symmetry one can compare e.g. with known values for critical exponents 
obtained by other methods [12, 46]. A comparison of our results for the critical exponents 
(3 and v with the results of the most sophisticated calculations show agreement within a 
few per cent. The anomalous dimension is also well determined, even though it is most 
affected by our truncation. 
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Finally, we should mention that our approach can be extended in several directions. 
The generalization to complex N x N matrices for arbitrary N is straightforward. For 
large N this opens the possibility of a comparison with 1/iV-expansions [47, 48]. Similarly, 
one may study systems with symmetry U{N, N) or U {2N) and symmetry breaking to the 
subgroup U{N) x U{N), relevant for the study of the metal insulator transition [1]. Very 
interesting generalizations are the systems with reduced SU (N) x SU (N) symmetry. They 
obtain by adding to the classical potential a term involving the invariant ^ = det (p+det ip^ . 
( Note that ^ is not invariant with respect to U{N) x U{N)). This will give an even richer 
pattern of phase transitions and permits a close contact to realistic meson models in QCD 
where the axial anomaly is incorporated. Finally one can extend the three dimensional 
treatment to a four dimensional study of field theories at nonvanishing temperature along 
the lines of ref. [6]. We hope to gain in this way new information about details of the 
chiral phase transition in QCD. 
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A Pole structure of the integrals 

The integrals 



I'nH = -nfdyy'^^'^[y{l + r{y))+u;]-^''^'\ 



^ ty^ POO ^ 



(A.l) 



with 

exhibit for d < 2(n + 1) a singularity at u; = —1. The massless dimensionless average 
propagator y{l + r(y)) is a monotonic function of y that takes on its minimum at 7/ = 
with limy^o2/(l + ^(y)) = 1- We define new variables 5 and 2;, 



5^u + l , z^y{l + r{y))-l (A.3) 



and substitute in (A.l), 



m = / dzGt{z)iz + S)-(-^'\ 
Jo 

POO 

li{8) = / dzGi{z){z + 5) 
Jo 



(n+l) 



(A.4) 



with 



ny2^ — — 



dy 



9y (A.5) 

ny2r{y) 



and y = y{z). For d < 2(n + 1) the integrals (A.4) have a pole at 5 = 0. (The singu- 
larity becomes logarithmic in 5 ior d — 2{n + 1)). In this case for 5 — > the dominant 
contribution to the integral comes from the region |/ ~ or equivalently z ~ 0. To find 
an approximate expression for and near the pole we expand the regular part of 
and Gi around z = 0. With 



r{y) = li^-ly + T^y' + oiy')), 



(A.6) 
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one obtains 



GHz) 

The inversion of z{y) given by (A. 3) can be done by expanding 



^yi-' (^2-^y + ^y' + 0(y') 



(A.7) 



and we find 



2.2 , 4 , 



(A.8) 
(A.9) 



Insertion of (A.9) in (A.7) yields 

GHz) = 2n(2z)i-^(l-i(^^ + l)z + 
GHz) = n{2z)'^-'(l-U^ + 4\z + 



1 + 

3 18 V 2 



+ 6 



z' + 0{z')j, 
z' + 0{z')] . 



(A.IO) 

We consider for d = 3 and n > 1 the zeroth order expression for and that obtains 
from the first term in (A.IO) and the exchange of summation and integration in (A. 4). 
Near the pole one finds 



roo 
^0 



-fni^)- 



(A.11) 



The leading contributions to If, /| and are therefore given by 

11(5) = 2V27rrV2, 

/3(<5) = 2-5/237r(5-5/2. 



(A. 12) 



We have verified this numerically. 
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